Resonant particle production in branonium 



OO 
O 

o 

(N 

< 

(N 



Or 

> 
oo 
in 

o 



o 



X 



J. G. Rosa^B and John March-Russell 1 '0 

1 The Rudolph Peierls Centre for Theoretical Physics, 
Department of Physics, University of Oxford, 
1 Keble Road, Oxford 0X1 3NP, UK 
(Dated: April 23, 2008) 

We study the mechanism of particle production in the world-volume of a probe D6-brane (or D6 
with SUSY breaking) moving in the background created by a fixed stack of D6-branes. We show 
that this may occur in a regime of parametric resonance when the probe's motion is non-relativistic 
and it moves at large distances from the source branes in low eccentricity orbits. This leads to an 
exponential growth of the particle number in the probe's world-volume and constitutes an effective 
mechanism for producing very massive particles. We also analyze the evolution of this system in an 
expanding universe and how this affects the development of the parametric resonance. We discuss 
the effects of transverse space compactification on the probe's motion, showing that it leads to the 
creation of angular momentum in a similar way to the Affleck-Dine mechanism for baryogenesis. 
Finally, we describe possible final states of the system and their potential relevance to cosmology. 
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I. INTRODUCTION 

One of the most important aspects of string the- 
ory with potential applications to beyond-the-standard- 
model physics and cosmology is the existence of addi- 
tional extended objects. A particularly interesting class 
of these, Z?p-branes, occur as (p+ l)-dimensional surfaces 
arising as supersymmetric solutions of the low-energy su- 
pergravity field equations but have also a high-energy in- 
terpretation as surfaces confining the endpoints of open 
strings [l| . By now many different configurations of Dp- 
branes have been studied, including parallel and inter- 
secting sets of branes 0, H|, the combination of branes 
with different dimensionalities [H and also several possi- 
ble embeddings in the compactified extra-dimensions Q . 
These models have an enormous potential in the con- 
struction of gauge theories with matter and several new 
extensions of the Standard Model have been obtained in 
this way 

Branes have also been intensively studied in their ap- 
plication to the problems of early universe cosmology, 
such as inflation and reheating. In the context of effective 
field theory treatments of brane-world model building, 
an immense amount of work has been performed on the 
astrophysical and cosmological implications of these the- 
ories for both flat and warped extra dimensions @ , while 
many studies in string theory have also been undertaken 
0)H)Q- One important property of Z?p-branes is the fact 
that they interact at large distances via the exchange of 
closed string modes, which allows one to study them as 
massive point-like charged objects moving along the di- 
mensions transverse to their world- volumes. An interest- 
ing application of this was proposed in (lfjj , where it was 



shown that branes with opposite charges may form bound 
states similar to electron-positron states, despite the dif- 
ferent origin of the central potentials involved. This anal- 
ogy has motivated the name of branonium to designate 
this type of systems. Several pro perties of these bound 
states have been analyzed in [lfj|, including the classi- 
cal orbits, the associated quantum dynamics and possi- 
ble compactification schemes. This analysis assumed a 
probe (anti)brane moving in the background spacetime 
created by a fixed central stack of branes, thus neglecting 
the gravitational backreaction of the probe. In the case 
of D6-branes, the probe's trajectory follows closed ellip- 
tical orbits similar to those found in planetary systems. 
It was also suggested that the probe's motion would lead 
to radiation of particles into both the bulk space and the 
probe's world-volume, making otherwise stable orbits de- 
cay. The stability of this system was further studied in 
[111 ] . where it was shown that orbital decay necessarily 
occurs when geometric moduli of the background space- 
time are allowed to vary. 

In this work, we explore some properties of branonium 
systems, namely the mechanism leading to particle pro- 
duction in the world-volume of the probe brane. For the 
particular case of a scalar field living in the world- volume 
of a probe _D6-brane, we show that particle production 
may occur in a regime of parametric resonance similar 
to the one behind the phenomenon of preheating after 
inflation [l2l [l3j . leading to an exponential growth of 
the particle number in the probe's world- volume.* We 
determine how the orbital parameters of the system con- 
strain this mechanism and analyze the possible proper- 
ties of the resonantly produced particles. Also, we study 
how the probe's trajectory is modified by compactifying 
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three of the dimensions parallel to the branes, so that a 4- 
dimensional spacetime with a flat Friedmann-Robertson- 
Walker geometry is obtained. In this case, we analyze the 
effects of Hubble expansion on the orbital parameters and 
how it modifies the parametric resonance regime. Other 
instability sources such as radiation into bulk modes are 
discussed as well. We also consider the effects of trans- 
verse space compactification on the probe's motion in the 
case of a compact 3-torus and analyze the associated pro- 
duction of angular momentum. Finally, we explore some 
of the cosmological applications of the resonant particle 
production mechanism and relate them to possible final 
configurations of the system, including stable states or 
brane-antibrane annihilation. 

In the next section, we begin by introducing some of 
the general properties of Dp— Dp branonia, with particu- 
lar emphasis on the p = 6 case, and discuss the effects of 
Hubble expansion on the orbital parameters. In Section 
III, we study the evolution of scalar particle-modes con- 
fined to the probe brane and analyze the properties of the 
associated parametric resonance mechanism in both the 
non-expanding and expanding universe cases. We dis- 
tinguish between the production of massless and massive 
particles and determine the significance of energy damp- 
ing into these modes on the probe's motion. We also 
consider, in this section, the effects of radiation into bulk 
modes on the resonant production of brane-bound par- 
ticles. In Section IV, we analyze how compactifying the 
directions transverse to the branes modifies the probe's 
trajectory and dynamically generates the required angu- 
lar momentum. Possible final states of the system are 
described in section V, where we also analyze the rele- 
vance of the particle production mechanism to cosmology 
in some particular scenarios. We summarize the main re- 
sults and conclusions of this work in section VI. 



II. EVOLUTION OF BRANONIUM IN AN 
EXPANDING UNIVERSE 

Branonium corresponds, in its simplest form proposed 
in jiq |. to a system whose components are a fixed stack 
of TV parallel p-branes and a probe p-brane, parallel to 
the stack, whose gravitational backreaction is neglected 
so that it has no influence on the geometry of the back- 
ground spacetime. Let the coordinates x^, fi = 0, • • • , p, 
correspond to the dimensions parallel to the stack and 
the coordinates y m , m = p + l, •• • , D —p, correspond to 
those transverse to the branes. Then, the general solution 
of the supergravity equations of motion in D dimensions 
for N parallel p-branes is given, in the Einstein frame, 
by [H: 



ds 2 



= h'^dxl 



f h^dy 2 
Cq\... p ■ 



ai-h- 1 ) 



(i) 



which give the D-dimensional metric, the dilaton $ and 
the non- vanishing components of the bulk (p + l)-form 



field Cfp+i] which couples to the (p + l)-dimensional 
world- volume of the stack of branes. The harmonic func- 
tion h(r) is a function of the radial coordinate in trans- 
verse space, r 2 = 6 mn y m y n , and has the form: 



h(r) 



Qp 



(2) 



where d = D — d — 2 and the factor of 1 ensures that the 
solution is asymptotically the D-dimensional Minkowski 
space. For the cases we will be interested in studying, the 
transverse space is at least 2-dimensional, so that d > 0. 
The stack of branes is, thus, located at the origin in the 
transverse space. The exponents 7, 7, k and the constant 
£ that characterize the solution are given by: 



Ad 



(D-2)A ' 

2a 

~A ' 



c = 



(D - 2)A 
2 



where d = p + 1 and 



2dd 
D-2 



(3) 



(4) 



The constant a n defines the coupling of the n — (p + 
2)-form field strength to the dilaton in the supergravity 
action. We will be interested in the particular case of 
Dp-branes in Type II supergravities, which arise as the 
low energy limit of Type II string theories, so that the 
?i-form field strength corresponds to a closed string state 
in the Ramond-Ramond sector, and a n = olr = 2(D — 
2n)/(D - 2). Taking D = 10, we have a R = ±(3 - p), 
A = 4, K = |(3 - p), d = 7 - p, 7 = |_(7 - p) and 
7 = \{p-\- 1)- It is clear that the condition d > implies 
p < 7. The constant Q p is related to the charge carried 
by the stack of N branes, and for the case of Dp-branes 
we have 



c P 9. 



,N11. 



(5) 



where c p = (2Y / 7r) 5_p r(^ £ ) , l s is the string length and 
g s is the string coupling constant. This solution holds in 
the limit of validity of string perturbation theory and for 
small curvature, which allows a low energy description of 
the system in terms of supergravity fields. This approxi- 
mation is valid as long as the local string coupling, g s e*, 
is everywhere small. Also, we need to require the radial 
coordinate r, which will give the distance between the 
probe and the stack branes, to be large compared to the 
string length l s , so that mediation by supergravity bulk 
fields (graviton, dilaton and RR-forms) is the dominant 
source of brane interactions. 

Before analyzing the motion of the probe brane, we 
need to consider possible compactifications of the 10 di- 
mensional spacetime that yield our 4-dimensional world 
at low energies. Compactification of dimensions paral- 
lel and transverse to the world volume of the branes 
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involve different procedures. In the case of dimensions 
parallel to the branes, it suffices to consider the Kaluza- 
Klein ansatz. If we start with the 15-dimensional space- 
time Eq. (TT]) and compactify along a periodic coordi- 
nate z ~ z + 2nR with radius R, so that the branes 
are wrapped around this dimension, we obtain a D' = 
(D — l)-dimensional spacetime with metric given by [l5j |: 



ds 2 D = e 2 ' aip ds 2 D , + e 2 ' bv {dz + B t ,dx f *) 2 
ds 2 D , = h-^'dx 2 p , + h 1 ' dy 2 , 



h" 



(6) 



where 



1 



2(D> - 1)(D' - 2) 

1 



b = -(£>' - 2)o 



p 2(d' - 2)a ' 
7' = 7 — 2p , 7' = 7 + 2ap . 



(7) 



The dilaton field $ is not altered by this compactifica- 
tion, while the (n — l)-form RR field can be decomposed 
as follows: 



C[n-i] = B[ n -x] + B[ n _2] A dz . 
We take the following truncation of these fields: 



(8) 



B 



[n-l] 



, 



(■B[n-2])oi 



ca-O- (9) 



This truncation leads to a (p — l)-brane supergravity 
solution in D — 1 dimensions, characterized by the same 
harmonic function of the transverse space radial coordi- 
nate h(r). Applying this toroidal compactification pro- 
cedure j times we obtain (p — j)-brane solutions living 
in a (D — j)-dimcnsional spacetime, with the number of 
transverse dimensions remaining the same. 

The process of compactifying the dimensions trans- 
verse to the branes is not as straightforward due to 
the lack of translational invariance which results from 
the presence of the branes. In particular, to make a 
given transverse coordinate z' periodic, so that a simi- 
lar Kaluza-Klein mechanism can be applied, one needs 
to include the appropriate "image branes" of the source 
branes which will ensure the invariance of the configura- 
tion under z' — > z' + 2irR' , R' being the radius of the 
compact dimension. For example, if a source brane is 
located at z' = 0, an image brane has to be considered at 
z' = 2t:R' . This implies a generalization of the harmonic 
function in Eq. © to include branes located at different 
points in transverse space [l5|: 



My) = i + E 



Qp 



(10) 



where the vectors yi denote the positions of the different 
image branes and |y - y ; | 2 = S mn (y m - yT)(y n ~V?)- 

These two compactification schemes allow us to reduce 
a 10-dimensional solution of the supergravity field equa- 
tions to a 4-dimensional spacetime which resembles our 



world. We will be interested in studying a particular case 
of the general branonium configuration described so far 
where we have a source stack of D6-branes and also a 
probe D6 brane. As was shown in [10|, when the probe 
carries an antibrane charge, i.e. opposite to the charge 
of the source branes, the system yields closed elliptical 
orbits which can be solved analytically, while for other 
values of p < 7 the orbits fail to close. This feature has, 
as we will describe later on, important consequences for 
particle production in the world- volume of the probe. 

Compactifying 3 of the dimensions parallel to the 
branes, we obtain effective D3-branes moving in a 3- 
dimensional transverse space, the shape of the orbits not 
being affected by this toroidal compactification scheme 
as mentioned before. If the transverse dimensions are 
finite but their typical size is much larger than that of 
the compact parallel dimensions and the interbrane dis- 
tance, we may consider them to be infinite and neglect 
the brane images described before. Although this would 
correspond to a 7-dimcnsional spacetime, from the point 
of view of fields confined to the world-volume of the 
source and probe branes it is effectively 4-dimensional. 
The dynamics in the transverse space will, nevertheless, 
affect the 4-dimensional dynamics of the fields. 

In order to make the system more similar to the observ- 
able universe, we consider a modification of the solution 
Eq. (1) so that, after compactification of 3 of the par- 
allel dimensions, the 4-dimensional world-volume of the 
branes has a flat Friedmann-Robertson- Walker (FRW) 
geometry. This is a particular case of the the system 
studied in [H| . where the moduli governing the size of 
both parallel and transverse dimensions were allowed to 
vary. In this work, we will assume that some dynamical 
mechanism at the string scale (see, e.g. [r|) or at some 
other high-energy scale (such as SUSY breaking) fixes all 
moduli except the scale factor associated with the 3 non- 
compact dimensions parallel to the branes. Hence, in the 
Einstein frame, the complete 7-dimensional line element 
is given by: 

ds 2 . = hT^l-dt 2 + a 2 6 l:j dx i dx 1 } + h'> 4 S mn dy m dy n , (11) 

where a = a(t), i,j = 1,2,3 and m,n = 1,2,3. The 
harmonic function is in this case h(r) = 1 + Qe/r. The 
exponents 74 and 74 can be obtained using Eq. ^ and 
applying Eq. |(7J| for the compactification on a 3-torus, 
but their precise values will not be necessary in the sub- 
sequent calculations. 

The motion of the probe brane is described by its ac- 
tion, which is divided in two parts. The Born-Infeld con- 
tribution, which concerns the geometry of the world vol- 
ume and possible gauge fields living on it, is in the string 
frame given by, for a generic p-brane probe [l(j | : 



~BI 



-T V J dP +1 i e~*^-det(% v + 2£?V) . ( 12 ) 



where £^ are world-volume coordinates, T p gives the p- 
brane tension and 7^^ = §MNd^x M d v x N is the induced 
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metric on the probe. The hatted metric tensor is defined 
in the string frame, and is related to the Einstein frame 
one via g^ v = e^g^ with A = 4/(D — 2) in the general 
case. The antisymmetric field strength tensor refers 
to gauge fields describing open string modes with end- 
points on the probe (in the case of Dirichlet branes), but 
in what follows we will set these fields to zero^ . The Wess- 
Zumino part of the action describes the minimal coupling 
of the Ramond-Ramond (p+ l)-form to the (p+ l)-world 
volume of the probe brane and can be written as [101 ] : 



probe brane's action to obtain, to lowest order, 



Swz 



-qT P / Ci 



y [p+i] 



(13) 



The constant q gives the sign of the charge of the probe 
brane, which being a BPS saturated state is equal or op- 
posite to the brane tension T p , so that q = 1 corresponds 
to a brane and q = — 1 corresponds to an antibrane [Iol |. 

If we now consider the particular case of D6 branes in 
10 dimensions, compactifying 3 of the parallel dimensions 
and including the scale factor describing the expansion 
(or contraction) of the 3 remaining non-compact dimen- 
sions along the brane, and choose the world-volume co- 
ordinates such that ^ = x^, we arrive at the following 
expression for the total action of the probe brane: 



T 6 V 3 / d 4 x a 3 (t) 



hv 2 



9 1 



(14) 

where V3 gives the volume of the compact 3-torus around 
which the D6-branes are wrapped to give the effective 
4-dimensional space. We have assumed that the trans- 
verse space coordinates are exclusively time-dependent, 
y m = y m (t) and defined the velocity of the probe as 
v 2 = S mn y m y n , with dots denoting time derivatives. As 
mentioned before, the bound orbits of the probe brane 
can be exactly solved for the case where the scale fac- 
tor is constant and it was shown in 10] that they cor- 
respond to closed ellipses. We now wish to analyze the 
effects of the varying scale factor on these orbits. In this 
case, a full analytical solution is hard to obtain and it 
is simpler to consider the evolution of the system when 
the probe brane moves at larges distances from the stack 
with small velocities. These two conditions can be ex- 
pressed as hv 2 <C 1. One must recall that we considered 
initially that the transverse dimensions can be considered 
infinite in extent. Hence, our analysis will hold as long 
as the radial coordinate is large compared to the typi- 
cal length scale of the harmonic function, given by Qe 
(assumed larger than the string scale and the compactifi- 
cation scale) , but still much smaller than the typical size 
of the transverse dimensions (assumed much larger than 
Qo). With this hierarchy in mind, we may expand the 



This means that we will consider the particular scenario where 
the classical background of these gauge fields vanishes. The role 
of gauge fields was considered in 



Sa 



T 6 V 3 / d 4 x a 3 (t) 



(15) 



Our subsequent analysis will not be affected by the 
constant term and we may drop it. For a probe antibrane, 



S 4 « T 6 V 3 I d A x a 3 



(*) 



(16) 



The rotational symmetry of the problem implies that 
the probe's trajectory will be confined to a plane, which 
we may choose to be y 3 = 0. Then, we may define the 
canonically normalized complex scalar field 



TqV 3 j . 2 



(17) 



so that its action is, in the non-relativistic and large dis- 
tance limit, 

S* = J d A x a 3 (t)( - g^d^dA* ~ V (</>)) , (18) 

where g^ v is the 4-dimensional flat FRW metric of the 
spacetime where the field is defined and the potential is 
given by: 



V{(j>) = 2T 6 V 3 h~ 



2T 6 V 3 [ 1 - 



nv 3 



2 |0| 



(19) 



so that it only depends on p = \(f>\. Note that in writing 
the action in this form one must take into account that 
the spatial derivatives of the field, di4>, i = 1,2,3, are 
assumed to vanish, i.e. the probe brane should remain 
parallel to the source branes during its motion. The anal- 
ysis of the classical motion of the probe brane is in this 
way reexpressed as a classical field theory problem of a 
complex scalar field in the background of a flat FRW 
spacetime. Before deriving the equations of motion for 
the field </>, it is useful to consider some of the symme- 
tries of the action. First, the rotational symmetry of the 
problem in the plane transverse to the branes is associ- 
ated with a conserved angular momentum. In terms of 
the field theory approach, this corresponds to the global 
U(l) invariance of the action under cj> — > e ia (f>, where a 
is a constant. The associated Noether current induces 
a conserved charge which gives, in the quantum theory, 
the particle number operator associated with the field. 
Writing the field as <fi — P e%e an d defining the comov- 
ing angular momentum as I = p 2 6, the comoving particle 
number density is given by: 



) = 2p 2 9 = 21 , 



(20) 



while the total particle number, which corresponds to 
the conserved Noether charge, is given by N = a 3 n. The 
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energy density and pressure of the field can be obtained 
from its energy-momentum tensor, yielding : 

e = T o = \<P\ 2 + V(4>) =p 2 + l - + V(p) , 

P 

V^ T ~^ = |0| 2 ^ Vtf) = (? + ^ - V(p) . (21) 
a z p z 



The covariant conservation of the energy-momentum 
tensor can then be expressed as dE = —pdV, where the 
energy of the field is defined as E = a 3 e and V = a 3 is 
the volume of the expanding (or contracting) flat FRW 
spacetime. Defining the Hubble parameter H = a/ a, we 
can write this in the form: 



P) 



(22) 



Prom Eqs. (JUl) and we conclude that ^ < 0, the 
equality holding only for H = 0. Thus, in an expanding 
universe (H > 0), the system evolves with a strictly de- 
creasing energy density. The angular momentum evolves 
according to: 



dl 
dt 



3Hl = Q 



(23) 



which corresponds to conservation of the total particle 
number AT. Hence, in an expanding universe, the ab- 
solute value of the angular momentum will necessarily 
decrease, vanishing asymptotically. 

We, thus, conclude from this simple analysis of conser- 
vation laws that the expansion of the universe will make 
the system reduce both its energy density and its angular 
momentum. To analyze the details of this evolution, wc 
need to determine the equations of motion for the field 
(j>. Varying the action Eq. (fl8|) with respect to <f>* , and 
taking this variation to vanish, we obtain: 



3H<j> 



8VM) 







(24) 



In terms of the transverse space coordinates y m , m 
1,2,3, this can be written as: 



y m + 3Hy m + 2Q 6 ^ T =0 



(25) 



This implies that the polar variables p and 9 satisfy: 

l 2 , \o_ 

p 3 ' 2 p 2 



6 + 2Z0 + 3H9 = 
where wc have defined the constant 



(26) 



equation for 9 simply gives the evolution of the angular 
momentum that we obtained previously in Eq. (|23[) . 

The evolution of the radial field p is determined by 
the form of the Hubble parameter, given by the energy 
density content of the universe via the Friedmann equa- 
tion in the usual way. This may include not only the 
energy density of the field <fi but also all other matter, 
radiation or vacuum energy components. The possibility 
that this scalar field may drive (slow-roll) inflation if, in 
the early universe, it dominates the energy density has 
been analyzed in [7|, for the case where the probe brane 
moves in a linear trajectory, i.e. I = 0. It was shown 
that the interbrane potential arising from Type II super- 
gravity /string theories is not flat enough to produce the 
required number of e- foldings, exhibiting an "^-problem" . 
This is intrinsically related to the assumption that r is 
much smaller than the size of the transverse dimensions. 

An inflationary period driven by the angular field vari- 
able 9 was proposed in where it is argued that, 
despite having a flat potential, this field may provide 
the constant energy density necessary for inflation if the 
probe brane moves in a circular orbit with p = 0. How- 
ever, we need to take into account the fact that, during 
inflation, the angular momentum decays exponentially as 
I oc e~ 3Ht , from Eq. (f2"3"| . The system will hence quickly 
tend to the I = case, where sufficient inflation is dif- 
ficult to obtain. Furthermore, the decay of the angular 
momentum will necessarily alter the evolution of p and 
the value of H, so that a slow-roll inflationary mechanism 
with the probe brane moving in a circular orbit would be 
hard to construct'''. 

Although we do not wish to completely discard such in- 
flationary mechanisms, we will from now on assume that 
some other field is responsible for inflation and analyze 
how the compactified D6 — D6 branonium system evolves 
in the post-inflationary eras. 

Before analyzing in detail the effects of the expansion 
on the trajectory of the probe, let us recall the general 
bound orbits of the system in a non-expanding universe. 
In this case, the system reduces to the well-known prob- 
lem of a particle in a central 1/p attractive potential, 
admitting closed orbit solutions of the form: 



P(0) 



R{l-e 2 ) 
1 + e cos 9 



(28) 



These are closed ellipses with eccentricity and semi-major 
axis given by: 



e = 1 



Ael 2 



2e 



(29) 



a = 4Q e 



.3/2 



(27) 



so that, apart from constant factors, the large distances 
potential is V(p) = ——. It is easy to check that the 



* Some work has been done recently in the context of B-branes 
moving in warped throats llgjl . In particular, the important role 
of the angular variables in providing accelerated periods of ex- 
pansion was discussed 



;s i n p 
in [3 
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The properties of such orbits may be obtained by 
studying the effective potential V e ff(p) = l 2 /p 2 + V(p). 
This potential has a minimum at — (a/21) 2 , tends to +00 
at the origin and to zero at large distances. The condi- 
tion p > then implies that orbits with —(a/21) 2 < e < 
will be bounded, with < e < 1. In this case, p 
will oscillate between its minimum and maximum val- 
ues, Pmin = i?(l — e) and p m ax = i?(l + e), which can be 
obtained from Eq. ([28]) by setting 8 = and 9 = tt, re- 
spectively. At the minimum of V e ff(p), the orbits will be 
circular with p = p c = 2l 2 /a, and e = 0. For e > 0, the 
orbits will only be bounded from below with e.g. p > I 2 /a 
for e = 0, which corresponds to an e = 1 parabolic orbit. 
The linear trajectory, with I = 0, is a particular case of 
the latter with no bounds on the value of p except for 
the trivial p > 0. For e > 0, the orbits are hyperbolic. 

Consider now the expanding case with the scale fac- 
tor evolving as a(t) oc t a . This power law behavior is 
typical of the post-inflationary stages of the universe's 
evolution, where one may consider a single fluid to give 
the dominant contribution to the total energy density. 
For example, in the radiation era we have a = 1/2 and 
in the matter era a = 2/3. We can also take a = 2/3 at 
the end of inflation, when the oscillations of the inflaton 
about the minimum of its potential, with an equation 
of state corresponding to that of non-relativistic mat- 
ter, dominate the energy density. This model only holds 
away from the transitions between these periods as in 
these cases at least two of the components give similar 
contributions to the energy density. 

We will consider the evolution of the branonium system 
for arbitrary a, starting at some instant to when the scale 
factor has a value an = a (to)- The Hubble parameter is 
then of the form H = j, which implies from Eq. |23|) 
that the angular momentum evolves according to: 



Substituting into Eq. (|3ip we obtain 



l = lo 



-3a 



(30) 



with l being its initial value. It is then clear that I — > 
as t — > +00, as we have mentioned earlier. The equation 
for the radial field p can now be written as: 



3a . 
P+ T P 



;2 

'-() 



-6a 



la_ 
2~p~ 2 



= 



(31) 



Let us start by analyzing how the system evolves when 
placed initially in a would-be circular orbit, so that 



2/ 2 



PcO 



(32) 



as p(to) = p(to) = 0. It is easy to see that p cannot 
remain constant for t > to due to the decay of the angular 
momentum. We may, however, look for solutions where 
the condition for circular orbits is maintained during the 
motion of the probe brane, i.e. 



Pit) 



2l 2 (t) 



- 6 a. 



PcO ( — ) = Pc(t) 



(33) 



PcO 
t 2 



-6a-2 



6a(3a + 1) = 



(34) 



which is satisfied only in the non-expanding case, a = 0, 
and for a = — 1/3. As we are interested in cases where 
a > 0, we conclude that p c (t) is not an exact solution of 
the equations of motion. It is, however, an approximate 
solution at late times and, as we will check later, gives the 
global evolution of the system in an expanding universe, 
so that it is useful to consider some of its properties. In 
particular, the angular frequency evolves as: 



I 



9o 



(35) 



so that the orbital period decreases as T c (t) cx (t/t )~ 9a - 
Such a solution eventually fails to satisfy the non- 
relativistic and large distance approximation as the ra- 
dius of the orbit decreases and its angular velocity grows. 
In particular, the parameter hv 2 which controls this ap- 
proximation grows, for t 3> to, as 



hv 2 



(36) 



Hence, given the values of T 6 and V3, it is the ratio Qe/lo 
that controls how long the approximations remain valid. 

Consider now trajectories with a non- vanishing initial 
eccentricity. In an expanding universe, the effective po- 
tential becomes time-dependent, although its asymptotic 
properties at the origin and at infinity remain the same. 
As the angular momentum redshifts away, the minimum 
of the effective potential decreases, the same happening 
with the value of p at this minimum, according to the cir- 
cular solution defined in Eq. ((33)) . As discussed earlier, 
the energy density of the field always decreases during 
the motion of the probe brane in an expanding universe. 
Hence, if the system is initially in one of the closed or- 
bits with e < 0, the latter condition will be satisfied for 
all t > to. This means that, at all times, the probe 
brane will be in one of the closed elliptical orbits defined 
earlier. From Eq. (|29f . we also conclude that the semi- 
major axis of the orbit will decrease as the energy density 
evolves to more negative values. The eccentricity of the 
orbits may, however, not remain constant, as it depends 
on the variation of both the energy density and the an- 
gular momentum. In fact, using Eqs. (f2"2"|) and l|2"3"|) , one 
can show that the eccentricity varies according to 

de 2 I 2 I 2 

— = -12H^r(3e + p) = -24H^r(2T + V). (37) 
at a z " eH 

where wc defined the kinetic energy of the field as T = 
p 2 + 4, so that e = T + V and p = T - V. 

Thus, in an expanding universe, for finite Iq, one ex- 
pects the eccentricity of the orbit to vary unless T = 
— V/2. This last condition is satisfied on average in the 
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non-expanding case, corresponding to the Virial Theo- 
rem. It does not necessarily hold, however, in an ex- 
panding universe, with the semi-major axis of the orbits 
decaying as discussed above. It is nevertheless clear that, 
at late times, as H and I decrease, the variation of the 
eccentricity of the orbits should be smaller. The sign of 
2T + V is also not definite, so that the eccentricity may 
either increase or decrease during the motion of the probe 
brane. 

From Eqs. and (J2S]), we obtain: 



R(t) 



Pc(t) 



1 



! (t) 



(38) 



Hence, the variation of the semi-major axis follows the 
decrease of the circular solution obtained earlier, being 
also affected by the variation of the eccentricity. 

In order to have a better understanding of the evolu- 
tion of the probe brane's motion in an expanding uni- 
verse, we have solved the equations of motion numer- 
ically. Measuring all quantities in terms of the string 
length, i.e. setting l s = 1, we choose, as an example, 
the values V 3 = (T 6 ) _1 = (2tt)- 6 and Q 6 = 100 § . Wc 
set the initial conditions at t = to to be those of a non- 
expanding elliptical orbit with angular momentum Iq and 
eccentricity eo, which we take to be the only free orbital 
parameters. The field is initially at its maximum value 



A) 



_ PcO 



j^f 2 -, with po = and 6 = n. After determining the 
numerical solution for p(t), we computed the correspond- 
ing energy density and eccentricity evolution, according 
to Eqs. (j2Tj) and ([29]) . respectively. We used the numeri- 
cal solution for the eccentricity to compute R(t), accord- 
ing to Eq. (|38p . and also the maximum and minimum 
values of p(t) at each orbit, given by: 



Pmin(t) = R{t)(l-e(t)) 

f J Tiiax 

(t) = R(t)(l+e(t)) 



(39) 



We have also determined the evolution of the angular 
field variable 9 using 9{t) = l(t)/p{t) 2 . 

In Figure 1 we have plotted the results obtained for 
the values Iq = 2000^ and e = 0.2 in a matter-dominated 
universe, a = 2/3, with t = 10 s **. 

Observing the plots shown in this figure, we conclude 
that, as expected, the probe brane evolves continually 
through elliptical orbits of decreasing semi-major axis, 
so that they fail to close. The radial field then oscillates 
between minimum and maximum values which decrease 
in time. The global decrease follows that of the circular 
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§ For the value of Qe we follow the example analyzed in 
■ The non-relativistic and large distance approximation holds for 
this value of initial angular momentum. 
** Notice that although this value of to seems quite large, one must 
take into account that it is measured in units of the string time, 
which is of order 10~ 43 seconds if the string scale is of the order 
of the Planck scale. Thus, we are considering the probe brane to 
be moving quite early in the history of the universe. 



FIG. 1: Numerical results obtained for lo = 2000 and e = 0.2 
in a matter-dominated universe, a = 2/3, with t = 10 s . The 
plots show (a) the radial field and associated quantities, (b) 
the energy density, (c) the eccentricity and (d) the motion of 
the probe brane in the complex plane of the field <f). 
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solution, Eq. ([33|) . but exhibits an oscillating behavior 
associated with the variation of the energy density. Al- 
though the latter strictly decreases, as expected, it also 
oscillates with decreasing period and amplitude, so that 
at later stages its evolution tends to be smooth. Conse- 
quently, one observes oscillations in the evolution of the 
eccentricity of the orbit. In the example shown in Figure 
1, these oscillations have a small amplitude. However, at 
earlier times, when the effects of Hubble expansion are 
more significant, numerical simulations show that this 
amplitude can be quite large. All simulations show that 
e(t) starts increasing and, as the period and amplitude 
of the oscillations decay, it tends to a constant value. 

A constant eccentricity at late times implies that, in 
this limit, the energy density should, on average, vary as 

( e (i))(xr 2 (t)(x (£) 6 " , (40) 

tending smoothly to more negative values as observed in 
Figure 1. Then, from Eq. (pj). 

(|) = -3H<e + p) = 6ff( e ), (41) 

so that {3e+p) = or, equivalently, (2T + V) = 0. 
Thus, although at early times the effect of the expansion 
makes the eccentricity vary significantly, at late times the 
system virializes and starts evolving smoothly between 
orbits of constant eccentricity and decaying radius. This 
result will be useful later for studying particle production 
in the probe brane at late times. 

Another interesting consequence of the expansion of 
the universe is a rotation of the axis of the elliptical orbits 
with time, although this is a small effect in the example 
shown in Figure 1 . This is simply a result of the growing 
eccentricity of the orbits which deviates the maximum 
and minimum values of p from 9 = n and 8 = 0, respec- 
tively. The orbital axis will, however, stabilize at late 
times, as e(t) tends to a constant value. 

Numerically, one also observes a decrease in the am- 
plitude of the oscillations of p{t). This is given by 

Ap(t) ee Pmax (t) ~ Pmin(t) = R(t)e(t) , (42) 

so that the observed decay of the amplitude is mainly 
due to the decrease of R(t) as discussed above. 

The asymptotic value of the eccentricity depends not 
only on the initial time to at which the probe starts its 
motion but also on the initial values of the orbital angu- 
lar momentum and eccentricity. Numerical simulations 
show that larger values of Iq lead to a more pronounced 
growth of the eccentricity, in agreement with Eq. (|37[) . 
The dependence on eo is less trivial and numerically one 
finds a larger eccentricity variation for eo close to and to 
0.9. This is, however, very small for eo close to 1, which 
is expected, as the energy density of the field cannot in- 
crease to produce hyperbolic orbits with e > 1 if initially 
e < 0. The eccentricity always grows if initially p — 0, so 



that 2T + V < and % > 0, as most of the eccentricity 
variation occurs initially when Hubble expansion is more 
significant. The opposite behavior should be observed if 
initially 2T + V >0. 

Although the example we have shown refers to a 
matter-dominated universe, a = 2/3, our qualitative dis- 
cussion and analytical results hold for all a > 0, in par- 
ticular to the a = 1/2 radiation-dominated universe. 

To summarize the results of this section, we conclude 
that, in an expanding universe, the probe brane follows 
elliptical orbits with decreasing radius and increasing 
frequency and whose eccentricity exhibits an oscillating 
behavior and asymptotically tends to a constant value 
larger than the initial one. If the motion of the probe 
brane begins at late times with a small angular momen- 
tum, the latter effect is negligible. 

The probe brane will asymptotically collide with the 
branes in the stack, as p — > 0, but our approximations will 
break down before this happens, so that we expect the 
motion of the probe brane to deviate significantly from 
our previous results as its angular velocity grows and it 
spirals towards the central stack. The collision will, how- 
ever, occur if nothing else prevents the probe from losing 
energy and angular momentum as the universe expands, 
and may lead to the annihilation of the probe antibrane 
with one of the source branes. For now, we will use the 
results obtained in this section to study the production 
of particles in the probe brane and we will return to this 
issue in section V of this work. 



III. PARTICLE PRODUCTION 

In the original discussion of the branonium system [lfj| | , 
it was argued that the fields confined to the probe brane, 
such as the gauge fields mentioned earlier, become time- 
dependent due to the motion of the probe through the 
background spacetime created by the central stack. Con- 
sequently, an observer bound to the probe interprets this 
variation as corresponding to the production of particles 
associated with these fields. Such radiation of energy 
into brane particle-modes arises only if the distance be- 
tween the probe and the stack varies in time, thus creat- 
ing a time-dependent background from the point of view 
of brane-bound observers. Although the power radiated 
into these modes was estimated in (loj . many aspects of 
this particle production mechanism remain unclear. In 
this section, we will analyze this mechanism in more de- 
tail, revealing some new properties of particle production 
in branonium, and discuss how the expansion of the uni- 
verse modifies this process. 

The fields confined to the brane arise, from the string 
theory point of view, as states associated with open 
strings whose endpoints are attached to the brane, as 
discussed earlier. For a single p-brane, which is a 1/2 
BPS supersymmetric state, these include bosonic fields 
such as gauge bosons and scalar fields, as well as their 
fermionic superpartners. If the probe corresponds itself 
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to a stack of M parallel p-branes, other degrees of free- 
dom associated with a U (M) vector supermultiplet may 
arise [20| . Here, we will study the simplest case of a real 
scalar field rj confined to the probe brane and follow an 
effective field theory approach. 

For the D6 — D6 system described earlier, we write the 
effective action for the scalar field rj in the form: 



S r] — —Tq 



J d 7 £ e-* 



where m denotes the mass of the scalar field, which in- 
cludes its bare mass as well as quantum corrections, and 
all other quantities are those defined earlier in this work. 
When compactifying 3 of the dimensions parallel to the 
brane as before, the scalar field will be decomposed in 
its Kaluza-Klein (KK) modes, leading to a tower of 4- 
dimensional massive modes, or KK states. We will focus 
on the evolution of the zero mode of the field rj, whose 
mass will simply be given by m. To simplify the notation, 
we will denote this mode as rj, although one must bear in 
mind that it is not the original 7-dimensional field. The 
global factor of Tq will not affect our discussion and we 
may absorb it into the definition of the field. 

First, let us consider the case m = 0. Then, the classi- 
cal equations of motion arising from the action Eq. (I43|) 
are given by: 



(44) 



Using the results obtained in Section II, and defining 
the function f(t) = [h(l — hv 2 )] 2, we can write this as 

fj + (3H - F)r) - ^(1 - hv 2 )V 2 rj = , (45) 



where F = j and V 2 rj = rj^didjT) is the flat 3- 
dimensional Laplacian of the field. 

To construct the associated quantum description, we 
follow the semi-classical approach to the quantization of 
scalar fields in curved space [Il[ . In a curved background 
spacetime with a time- varying geometry, as is the case of 
the world-volume of the probe brane, the induced time- 
variation of the fields modifies the usual canonical quan- 
tization procedure, as the quantum operators become 
themselves time-dependent. In particular, the creation 
and annihilation operators associated with the field will 
now evolve in time as the background changes, the same 
happening with the associated multi-particle states. This 
is the main reason behind the production of particles in 
a dynamical background. Let us start by expanding the 
field in Fourier modes of the form: 



77(x, t) 



d 3 k 
(2^)1 



(a k x k (t)e ikx + 4x£(i) e - ikx ) . (46) 



In the quantum theory, ak and a k become the annihi- 
lation and creation operators associated with the Fourier 
mode k of the field. Expanding the field in this way, we 



include all the time dependence of the field in the func- 
tions Xk(i) and X k (*)i while the operators remain time- 
independent. The conjugate momentum to the field rj, 
obtained by computing SS v /5r), can then be written as 



tt(x, t) 



d 3 k 

f J (^)§ 



/ • zk-x 1 t ** — tk-x\ / ai-t\ 

(akXke +a k Xk e )■ (47) 



If the creation and annihilation operators satisfy the 
canonical commutation relations 

[ak, a v ] = [al 4,] = , [a k) a£,] = 6 3 {k - k') , (48) 

then the canonical commutators for the field and its con- 
jugate momentum, [r?(x, t), 7r(y, t)] = ii5 3 (x— y), can only 
be obtained if the mode functions satisfy the following 
Wronskian normalization condition: 



(49) 



This condition is essential for consistency of the quan- 
tization procedure. We now follow a semi-classical ap- 
proach and consider the evolution of the field modes to 
be given by the classical equations of motion, assuming 
that any quantum corrections to their propagation can 
be neglected. Substituting Eq. (|4*6"j) into Eq. (|4"B")) , we 
see that each of the Fourier modes evolves independently 
according to 



Xk + (3H - F) Xk + —(1 - hv 2 )x^ = 



(50) 



and similarly for the complex conjugate modes xL- No- 
tice that this equation depends only on k = |k| but not 
the direction of the momenta k = ¥• This is a conse- 

k 

quence of the isotropy of our metric ansatz, and we may 
write without loss of generality Xk(t) = Xk(t)- It is also 
clear, from Eq. ([5U)) . that the physical momentum of the 
modes is related to their comoving momentum k via 



^phys 



'-Vl-hv 2 



(51) 



This expression includes both the expected redshift of 
the modes in an expanding universe and the effect of the 
probe's motion through the dynamical background. 

Let us now rescale the field modes by defining the mode 
functions Xk(t) = a? (i)/~ 5 (t)xk(t)- These satisfy 



X k + Lo 2 X k = 



(52) 



which corresponds to the equation for a harmonic oscil- 
lator with a variable frequency given by: 



<4 = 1 ^ I (l~hv 2 )-\{ZH-F) 2 - l -{3H-F) 
- h 2 4- A 2 

— K phys + ^ ) 

where we have defined 
1 



(53) 



^m-FY--{m-F) . (54) 
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As each mode behaves as a harmonic oscillator, we can 
define the associated particle number via 



E, 



n k 



X k \ 2 + -co 2 k \X k \ 2 



2 



\X k \ 2 



uj k [ n ! 

1 

2 



1 



(55) 



It is well known that this quantity is an adiabatic in- 
variant for a harmonic oscillator with variable frequency, 
as happens, for example, in an oscillating pendulum 
whose length is decreased infinitely slowly [2l|. Thus, 
the number of quanta in a given Fourier mode k can 
only change if its frequency varies in a non-adiabatic way, 
which can be expressed by the following condition [l3| : 



du k 



dt 



(56) 



For our particular case, this gives, in the limit hv 2 <C 1, 



\Hk 2 phys - AA| > (k 2 



phys 



(57) 



This condition will then constrain the number of parti- 
cles produced in each physical momentum mode by the 
motion of the probe brane. 

To quantify this number of particles, we need to look 
for solutions of the mode equation, Eq. |52|) . We will 
start by discussing what happens in a non-expanding 
universe and then we will analyze the effects of Hubble 
expansion on particle production. 



A. Non-expanding universe 

Consider the motion of the probe brane in a non- 
expanding universe, H = 0, and assume that the non- 
relativistic and large distance approximations discussed 
earlier are valid. As we have seen in Section II, the probe 
follows closed elliptical trajectories if its energy (density) 
is negative, e < 0. If hv 2 < 1, we have to leading order: 

1 - hv 2 ~ 1 , k phys ~ k , f~ h 1 / 2 . (58) 

Then, we may take the following approximations: 



F 



F 



so that we can write: 



-F 2 



f 
f 
1 
2 



-F 



lh 
2h ' 
h 
h 



5 

16 



1 h 

Ah 



(59) 



(60) 



An exact analytical expression for the time variation of 
the radial field pit) can be obtained in a non-expanding 
universe, being given implicitly by 



p(t) =R[l- ecos (V>(*))] 



(61) 



The angular variable ip(t) satisfies 

t — to = S! -1 Up — ?po + e(sin ip — sin?/>o)] 



(62) 



where ipQ = ip(to) and Q — \/jR~^ is the angular fre- 
quency of the orbit. The angular variables ip and 6 are 
related by: 



,'1 + e ^ 

tan — = \ tan — 

2 V 1 - e 2 



(63) 



For small eccentricity orbits, e« 1, Eq. (|62p gives, to 
leading order in e: 

V> ~ + n(i - t ) • (64) 

Choosing ip$ = 7r we have, in this approximation, 

p(t) ~ R[l + ecos (Q(t - t Q ))] . (65) 

This allows us to determine the non-relativistic and 
large distance expression for A 2 , in the case of small ec- 
centricity orbits. To lowest order, we then obtain: 



A^ 



1 T 6 V 3 Q ( 



en 2 cos (fi(t - t )) 



Defining 



nv 3 q ( 



R 



(66) 



(67) 



the variable frequency of the mode with comoving mo- 
mentum k can, hence, be approximated by: 



wl(t) ~k 2 +S 2 cos (n(t-t )) 



(68) 



If we now rewrite Eq. (|52|) in terms of the rescaled 



time variable z = U(t — to), we obtain: 



X k / +(A k -2qcos(2z))X k = Q, 



where 



A k 



S 2 



(69) 



(70) 



E q. ([551 has the form of the well-known Mathieu equa- 
tion [22j with parameters A k and q determined by the or- 
bital parameters of the probe brane and by the value of 
the comoving momentum k of each mode. In the param- 
eter space (A k ,q), the Mathieu equation exhibits both 
stable and unstable solutions, the latter being closely as- 
sociated with the phenomenon of parametric resonance 
[23l ] . As we are working under the assumption that the 
probe brane's orbit has a large radius and a small eccen- 
tricity, it is easy to check that \q\ 1 (the analysis is 
independent of the sign of q) . In this region of parameter 
space, the Mathieu equation exhibits instabilities in a se- 
ries of narrow resonance bands near A k ~ n, with n being 
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a positive integer, and with a width in comoving momen- 
tum space approximately given by Ak n ~ \q\ n - In these 

resonance bands, the solution evolves as X k oc e M * >z , 

(n) 

with a real exponent fi k , according to Floquet's Theo- 
rem. 

The most important of these is the first resonance 
band, which occurs when 



1 



1 



l-M-779 <A k <l + \q\--q< 



8 



(71) 



The exponent fik = fJLL can then be approximately 
written as 



~vV-04 fc -l) 5 



(72) 



This is, as claimed, real for 1 — \q\ <Ak < 1 + hav- 
ing a maximum value of \q\/2 for Afe ~ 1. This implies 
that particle-modes lying inside this resonance band will 
be exponentially amplified, leading to a resonant pro- 
duction of particles. Thus, from Eq. ([70]) . the center 
of the resonance band occurs for comoving momentum 
k c ~ f2/2, while its lower and upper bounds are given by 




6 -^k c 
2 



1- ^ 



1+ ^ 



(73) 



so that the resonance band has a width in comoving mo- 



mentum space of Ak res = k Ti 



krt 



This is 



quite small compared to the value of the center of the 
: \q\ <C 1, which means that parti- 



Ak r 



resonance, as 

cle production occurs in a regime of narrow parametric 
resonance. The particles produced by this mechanism 
have typical energies of order ui k — §, i.e. with half 
the typical energy of the field <fi. This agrees with the 
discussion in where it is stated that particles are 
produced in pairs with opposite momenta, as if resulting 
directly from the decay of the interbrane distance field' 
We have found in this work, however, that these parti- 
cles are produced resonantly, the associated field modes 
being exponentially amplified. 

After the modes inside the resonance band have been 
amplified sufficiently, we may write, from Eq. (|55p . 



n k cx e 2 ^ kZ = e »* n (t-t ) 



where we used that (if- < ^ <C 1, so that the term in- 
volving \Xt\ 2 in Eq. (f5"5")) can be neglected. This also 
implies that, although the particle number in each mode 



(74) 



tt This field can be seen as a classical condensate of zero momen- 
tum particles, so that only pair production of r?-particles ensures 
momentum conservation. 



grows exponentially, the resonance takes a long time 
to develop. In particular, the typical resonance time, 
At res = (/ifcfi) -1 , is much larger than the orbital pe- 
riod of the probe brane, which is of order 27rf2~ 1 . This is 
characteristic of the narrow resonance regime and, hence, 
a significant particle number is only produced after the 
probe brane has completed a large number of orbits, 
given by (2irfik)~ 1 - 

To have an idea of the order of magnitude of these 
quantities, consider the scenario studied in Section II, 
where V 3 = Tq 1 = (2tt) 6 with Q 6 = 100. Then, for 
an orbit with eccentricity e = 0.01 <C 1 and an angu- 
lar momentum I = 500, we have ^/T 6 V" 3 /2% ~ 0.03, 
which satisfies the large distances approximation, and 
\q\ ~ 10~ 4 <C 1. This implies that the probe brane will 
have to complete approximately 3 x 10 3 orbits for the 
particle number to be amplified by one e-folding. 

To illustrate the results obtained earlier, we have 
solved Eq. ([52")) for the example described above. An ini- 
tial vacuum state with zero particle number corresponds 
to the following initial conditions: 



X k (t ) 



1 



i/2wfc(io) 



X k (to) 



UJk(t ) 



(75) 



In Figure 2 we have plotted the numerical results ob- 
tained for the comoving momentum at the center of the 
resonance band, k c . The linear evolution observed for the 
logarithm of the particle number shows that, as expected, 
this quantity is being amplified exponentially with time, 
but that a significant number of particles is produced 
only after the probe brane has completed a few thou- 
sands of orbits. 




FIG. 2: Numerical results obtained for the evolution of the 
particle number associated with the mode in the center of 
the resonance band for an orbit with e = 0.01 and I = 500, 
giving q ~ 10~ 4 . The time coordinate t is given in units of 
the orbital period T = 2 7 rfi~ 1 ~ 1.57 x 10 5 . 



We note that, although this solution was obtained nu- 
merically, a formal analytical solution of the Mathieu 
equation can be written in terms of the Mathieu sine 
and cosine functions. 

In Figure 3 we also show the evolution of the particle 
number for a mode outside the limits of the resonance 
band. In this case, [i k takes imaginary values and n k 
exhibits an oscillating behavior with a small amplitude, 
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confirming that only for comoving momenta k m i n < k < 
k-max a significant particle number is produced by the 
elliptical motion of the probe brane. It is also worth 
mentioning that, in the limit of circular orbits, e — * 0, 
the parameter q vanishes and there is no amplification 
of any of the particle momentum modes. This confirms 
our early assumption that the interbrane brane distance 
needs to vary for particle production to occur. 




10000 15000 20000 25000 30000 



FIG. 3: Evolution of the particle number for a particle-mode 
outside the resonance band for an orbit with e = 0.01 and 
I = 500, giving q ~ 10 -4 . The comoving momentum of the 
mode is k — k c (l + 4^-) > k max . The time coordinate t is 
given in units of the orbital period T = 2irQ.~ 1 ~ 1.57 x 10 5 . 



So far we have not included the backreaction effects of 
the produced particles on the probe's motion. We will 
return to this later in section III C, where we estimate 
the energy density damped into brane-modes by the res- 
onance and discuss how this affects the probe brane's 
orbital decay. 

Consider now the case where the field has a non- 
vanishing mass, m 7^ 0. In this case, the classical equa- 
tions of motion Eq. (|44|) can be generalized to 



or, explicitly, 

77 + (3H - F)fj - -g (1 - hv 2 )\7 2 v + 



7m 2 77 







(1 - to 2 ) 



(76) 



m rj = 



(77) 

Each of the Fourier modes then evolves independently 
according to 



X k + (3H-F)xk- 



k 2 

{l-hv 2 )% + 



(1 - hv 2 ) 



(78) 



This implies that all modes have, in this case, an ef- 
fective physical mass which is different from the mass 
parameter m due to the motion of the probe brane: 



m lff 



(1 - hv 2 ) 



(79) 



Notice that the mass is not, however, affected by the 
expansion of the universe. When writing the equations 
of motion in terms of the rescaled mode function (t) 



defined earlier, we again obtain the equation for a har- 
monic oscillator with a varying frequency, which is now 
given by: 



= k 



phys 



A 



(80) 



where all the quantities are time-dependent in the gen- 
eral case. Let us now set H = and compare these 
results with the ones obtained for the particle produc- 
tion mechanism in the massless case. First, notice that, 
in the non-relativistic and large distance approximation, 
we have ni 2 ^ ~ m 2 . Thus, the leading order modifica- 
tion introduced by the mass of the field can be obtained 
by simply replacing k 2 — > k 2 + m 2 . 

Following the same procedure as before, one can reduce 
the particle-mode equation to the Mathieu equation, with 
the parameter Ak being now given by: 



4,(k 2 + m 2 ) 

n 2 



(81) 



The parameter q, which quantifies the strength and 
width of the resonance phenomenon, is not modified by 
the introduction of the mass parameter, so that we expect 
the particle production mechanism to occur as before. 
In this case, however, the center of the resonance band 
is shifted to a lower momentum, k c = <J (Vl/2) 2 — m 2 , 
which gives Ak = 1. The upper and lower limits of the 
resonance band are also modified, being now given by: 



kn 

kn 



(Q/2) 2 -m 2 ~ i<5 2 , 
{Q/2) 2 -m 2 + ^5 2 , 



(82) 



which in turn alters the width of the resonance band 
in the obvious way. If, as in the massless case, we have 
S 2 /2 <C k 2 , we may write Ak res ~ S 2 /2k c , which is larger 
than the corresponding value in the massless case, satis- 
fying nevertheless the narrow resonance condition. No- 
tice that this implies an upper bound for the mass of the 
particles which can be produced by this mechanism, 



m < 



'ny _ 1 2 _ n 
.2 2 ~ 2 



(83) 



In the particle description discussed earlier, this simply 
means that each "particle" in the interbrane distance 
field 4> needs to have enough energy to decay into two 
77-particles at rest. A residual kinetic energy is also nec- 
essary if all the modes inside the resonance band are to 
be excited. In the non-relativistic and large distance ap- 
proximation, the orbital frequency of the probe brane is 
assumed to be small, so that only particles with small 
masses may be excited during its motion. One must 
recall, however, that we are measuring all quantities in 
terms of the typical string parameters. In particular, the 
orbital frequency should be small compared to the string 
energy scale. If the latter corresponds to the Planck scale, 
then the produced particles can still be quite massive. 
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For the example studied earlier in this section, with 
orbital parameters J]~4x 10 -5 , particles up to masses 
of order 10 13 GeV can be resonantly created. Recalling 
the definition of <f> in terms of the physical interbrane 
distance, Eq. (|17[) . we may write the upper bound on 
the mass of the produced particles in terms of the phys- 
ical semi-major axis of the orbit, R p hy S — y^/ToVsR. 
Inserting back the missing l s factors, we obtain: 



\F\ > 3H and \F\ > |3i?|, or explicitly 



m,, 



Qe 



I, 



2R 



phys Rphys 



(84) 



where M s = l^ 1 is the string energy scale. As we assumed 
in the beginning that r 3> Qq 3> l s , it is clear from this 
expression that m max <C M s , which allows nevertheless 
for the production of very massive particles. In general, 
particles living in the world-volume of the probe brane 
will be massive, so that we need to use the results ob- 
tained above to describe the parametric resonance. For 
masses which are small compared to the typical orbital 
energies of the probe, the simpler treatment of the m = 
case will be sufficient. 



B. Expanding universe 



Let us now consider the case of an expanding uni- 
verse, with H = j > 0. The inclusion of the scale 
factor will modify our previous analysis of particle pro- 
duction in several ways. First, it will change the fre- 
quency of each harmonic oscillator particle-mode, namely 
modifying the factor A 2 (i), given in Eq. ([54")) . Also, 
the physical momenta of the particle-modes will be red- 
shifted as the universe expands, so that, in the non- 
relativistic and large distance approximation, we have 
kphys — k/a oc (t/to)~ a . Finally, the motion of the probe 
brane will be altered, as we have seen in Section II, so 
that the semi-major axis of elliptical orbits will decrease 
as R(t) cx (t/t )~ 6a , making the angular frequency con- 
sequently increase as D,(t) oc (t/t Q ) 9a . Furthermore, if 
the Hubble parameter is sufficiently large, the orbital ec- 
centricity may exhibit a significant growth. 



1. Analytical results 

All the modifications make a complete analytical study 
of particle production in the probe brane rather complex. 
In order to determine the leading effects of the expan- 
sion, we will consider the motion of the probe brane at 
sufficiently late times, so that some of the effects of the 
expansion can be discarded. As we have concluded in 
Section II, the variation of the eccentricity of the orbits 
is, in this regime, very small, so that it may be neglected. 
The explicit effects of the expansion in modifying the fac- 
tor A 2 may also be discarded in this limit, provided that 



6a f R phys \ 



t 2 > 



6a 



R 



phys 
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The first of these conditions will, a priori, be more con- 
straining. For the example we have been following in this 
section, the first condition gives t » 5 x 10 8 while the 
second only implies the constraint f > 3.5 x 10 6 . 

With these constraints in mind, we now have to analyze 
how Hubble expansion modifies the factors F and F in 
the variable harmonic oscillator frequency. The radial 
interbrane distance field p can now be written, in the low 
eccentricity orbit approximation, as: 



p(t)~R(t)[l 



(n(i)(*-t ))] 



with 



R(t) = R [- 



-6a 



n(t) =fl ( - 



(86) 



(87) 



Computing its time derivatives, we find that they in- 
volve terms arising from Hubble expansion which are sup- 
pressed by at least one power of t, being negligible at 
late times. There are also terms involving the quantity 
. If the time necessary for particle production 



At 

to 



_ t-t 

to 



to occur is smaller than the age of the universe at that 
time, we may also discard these terms. After some alge- 
bra, the leading order modification to A 2 is then given 
by 



A 2 (t) ~ S 2 (t) cos(fi(t)(t-t )) 



The quantity S 2 is now time-dependent and can be ob- 
tained trivially by including the appropriate time varia- 
tion of the orbital parameters in its non-expanding ex- 
pression, Eq. (|67|) . This gives 



1 T,V„Q, 



2 flo 



(89) 



Then, taking into account the momentum redshift, we 
can write the particle-mode equation, to lowest order, as: 



-la 



24a 



cos (n(t)(t-t )) x k = o. 

(90) 

Following the same reasoning as in the non-expanding 
case, we rescale the time coordinate by defining the vari- 
able z = ^ip-(t — to). As the angular frequency of the 
orbit is now time-dependent, this change of variables will 
introduce new terms in the equation. These will, how- 
ever, be suppressed by powers of t and involve first and 
second derivatives of the mode function X k ■ In the non- 
expanding case, these were quite small compared to X k , 
and one expects them to be negligible in this case as well. 
Using that, to lowest order, z ~ %Ai, we may write the 
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leading order equation for the mode function in terms of 
the variable z as 



X'l + [A k (z) - 2q(z) cos(2z)]X k = , (91) 



with 



A k (z) = A k0 (l-~fz) 



q(z) = qo(l+&), (92) 



where A k o = 4fc 2 /Oq and qo = ~25 2 /fl 2 1 correspond to 
the initial values of the parameters of the Mathieu equa- 
tion, given by their non-expanding case expressions, and 
the coefficients 7 and £ are given by 



7 



40a 



12a 



To 7 



(93) 



We, hence, conclude that the leading order effect of 
the expansion is to make the coefficients of the Mathieu 
equation vary in time. In particular, their variation is, 
to a first approximation, linear in the variable z and con- 
trolled by the single parameter 7, which will be small at 
late times. One must take into account, however, that 
Eq. (|9"Tj) is not the Mathieu equation and, namely, does 
not satisfy the conditions of Floquet's Theorem which 
allowed us to define the coefficient \x k . Nevertheless, we 
can use the results obtained in the non-expanding case 
from the properties of the Mathieu equation as a guide 
to describe the behavior of the modes when H > 0. 

First, notice that, as A k decreases with z, the position 
of the center of the resonance band, at A k — 1, will be 
shifted to higher momentum modes: 



k c (z) 



V<;0 



1 



7Z 



7 _ ^0 

kco = — 



(94) 



Thus, we expect modes with comoving momentum larger 
than fc c o to be excited as the probe's orbit decays. To bet- 
ter understand the evolution of each mode, let us analyze 
the evolution of the exponent fi k (z), which is given by: 

/**(*) - 5 V?o 2 (l + - (>M1 - 7*) ~ l) 2 • (95) 



The zeros of jJL k (z) occur for 



z\ = 



-jgbj + A kQ - 1 
7-^fco + £|<?o| 



_ M + Afco-i 
Z2 ~ — a c\ — r ' 

lA m ~ g\q Q \ 



with Z2 > z\. This means that [i k (z) will be real for 
Z\ < z < z 2 and pure imaginary for z < Z\ and z > z 2 . 
Figure 4 illustrates the typical evolution of /i k (z) for a 
generic Fourier mode. 

Thus, each mode may only experience the resonant 
regime during a finite amount of time, when /J, k (z) is real 
and the mode is inside the resonance band. During the 
periods where /i k (z) is pure imaginary, we expect the 
particle number to exhibit an oscillating behavior. In 
the non-expanding case, the frequency and amplitude of 
these oscillations for modes outside the resonance band 
is determined by the value of Im ji k , so that the larger 



Im iu k \ Re /u k / 




FIG. 4: Evolution of the real and imaginary parts of the expo- 
nent Hk(z) for a generic mode. The shaded area corresponds 
to the period the mode spends inside the resonance band. 



the latter the smaller the period and the amplitude of 
the oscillations. As, in an expanding universe, Im n k (z) 
changes in time, we expect the frequency and amplitude 
of the oscillations of the particle number to vary as the 
probe brane's orbit decays. Similarly, as Re fi k is not 
constant when the mode is inside the resonance band, 
the strength of the resonance is also expected to vary. 

Particle modes will then have different behaviours ac- 
cording to their position relative to the initial resonance 
band, which can be parametrized by: 



A w = 1 + /% | • 



(97) 



Recalling that the centre of the original resonance band 
corresponds to (3 = and that the probe's motion begins 
at z = 0, we need to distinguish three different types of 
modes: 

I. p < -1 

As both z\ and z-i are negative in this case, these 
modes will always be outside the resonance band, 
with an imaginary exponent fj, k (z) for all z > 0. 

II. -1 < (3 < 1 

This case corresponds to the particle-modes inside 
the resonance band in a non-expanding universe, 
with z\ < but Z2 > 0. Hence, they will start 
inside the resonance band but stop experiencing the 
resonant regime for z > 

III. /3 > 1 

Although these modes are outside the initial reso- 
nance band, they will experience the resonant be- 
havior at some later stage, as z 2 > z\ > 0. 



For both type II and type III modes, the maximum 
lue of 1 
8b I < 1, 



value of Re /i k (z) occurs for z max = Zl ^ Z2 , giving, for 



VkyZmax) ^ — I 1 + Y^P|<?0 



(98) 
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This gives the same result as in the non-expanding 
limit for (3 — 0. It is easy to check that z max -» ±00 as 
7 — > 0, the same happening for z\ and z 2 . This implies 
that only modes for which z\ < (— > —00) and Z2 > 
(— > +00) will be excited in this limit, in agreement with 
the fact that only type II modes are excited in the non- 
expanding case (7 = 0) and that these experience the 
parametric resonance during an infinite period of time. It 
is also clear from Eq. ((98)) that higher momentum modes 
are excited by the resonance at later times, as expected 
from the evolution of the center of the resonance band. 

The total time each type II or type III mode spends 
inside the resonance band is given by: 



Az 6 , 



and 



Z2 - Z\ 



gjgbj 

7 



1 



(99) 



It is clear that this quantity decreases with 7, so that 
the smaller the Hubble parameter the more time each 
mode is excited by the resonance. In the limit 7 — > 0, 
^ z band — * +00, but as previously discussed this limit 
only applies to type II modes, as all the others are outside 
the resonance band in this limit. One also concludes that 
higher momentum modes will spend less time inside the 
resonance band, which suggests these will be less excited. 

The analysis of the properties of the exponent /ik{z) 
thus provides a very clear insight on the qualitative evo- 
lution of all momentum modes. We would like, however, 
to be able to quantify the particle number produced by 
the resonance in each mode. In the non-expanding case, 
Xk(z) oc e MfcZ , which implied n^(z) oc e 2MfeZ in the limit of 
large particle number for type II modes. Despite the vari- 
ation of fik in the H > case, we may take this exponent 
to be approximately constant during an infinitesimal in- 
terval between z and z + dz. Then, during this interval, 
Xfc will be amplified by a factor exp(^ik(z)dz) if the mode 
is inside the resonance band. Integrating this result, we 
expect that, after a significant number of particles has 
been produced, 



n k (z) oc exp ( 2 / fi k (z')dz' 



(100) 



The initial time Zi refers to the time the mode enters 
the resonance band, i.e. z% = and z, = z\ for type II 
and type III modes, respectively. One can use Eq. (|100[) 
to estimate the total particle number produced in each 
mode by the parametric resonance. The details of this 
calculation are given in Appendix A. For type II and type 
III modes, one obtains: 



logn fc >{z 2 ) 



27l+/%o| 
arcsin/3) , 

27l + /%o| 



(101) 



Eq. (|10ip constitutes the main result of this section, 
giving the leading order expressions for the particle num- 



ber density produced by the resonance in an expand- 
ing universe. From these approximated expressions, we 
can see that the parameter Qg/7 controls the strength 
of the resonance, a significant number of particles be- 
ing produced only if 7 < qfi. Recall that, in the non- 
expanding case, the typical resonance time was given by 
At res — (/ifcS!) -1 ~ (|go|^o/2) _1 - The typical time a 
mode spends inside the resonance band in the expand- 
ing case is, from Eq. (|99)l . of order At}, an d ~ -^^q 1 . 
Hence, apart from numerical factors, the condition for a 
significant number of particles to be produced is simply 
stating that the time a given mode spends inside the res- 
onance band should be greater than the typical time the 
resonance takes to develop, Atband > Ai res . 

From the results in Eq. IjlOip . we can also conclude 
that the total particle number produced increases with 
(3 for -1 < /3 < 1. On the other hand, for (3 > 1, it 
exhibits a slow decrease with (3, showing that the pro- 
duction of high momentum modes is suppressed. This 
agrees with the fact that high momentum modes spend 
less time inside the resonance band, as obtained above. 



2. Numerical simulations in an expanding universe 

To check the results obtained so far for the resonant 
particle production in an expanding universe, we have 
solved Eq. (f9"Tj) numerically. In order to simplify the 
computation, we have rescaled the mode functions via 
Afc = y/2LJk(to)Xk, so that initially we have approxi- 
mately 



X k (0) ~ 1 



**(0) 



We then computed the function 

n fe (z) = |A fc (z)| 2 + |Al(z)| 2 



(102) 



(103) 



This is related to the physical particle number density of 
each mode, given by Eq. (|55]l . approximately by n& — 
, taking the limit of large particle number. Note 

that fit = 2 initially, as we neglected the 1/2 factor in 



Eq. (|55j) . so that the effects of this normalization have 
to be taken into account before the resonance produces 
a significant number of particles. 

In Figure 5, we illustrate the results obtained for modes 
of types I, II and III, including also the analytical pre- 
diction obtained from Eq. (|100|) . These simulations cor- 
respond to an original resonance strength |qo| = 10~ 3 , 
the effects of the expansion being quantified by 7 = 10~ 7 
(recall Eq. (f9"3"|) ). These values ensure the validity of the 
approximations made and also that a significant particle 
number is produced for modes of types II and III. 

Observing the plots shown in this figure, we see that all 
three particle-modes follow the expected behavior. The 
type I mode exhibits an oscillating particle number, with 
oscillations of decreasing amplitude and period, confirm- 
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FIG. 5: Numerical results obtained for the particle number 
with \qo \ = 10~ 3 and 7 = 10~ 7 . The plots correspond to (a) a 
type I mode, with /3 = —3, (b) a type II mode, with /3 = and 
(c) a type III mode with f3 = 3. The solid line corresponds to 
the numerical solution in all three cases, while in (b) and (c) 
the dashed line gives the corresponding analytical prediction. 
The shaded area corresponds to the period that each mode 
spends inside the resonance band. 



ing that, in fact, no net particle number is produced^. 
Both the type II and the type III modes are exponentially 
amplified during a finite period, after which their particle 
number oscillates with decreasing amplitude and period, 
tending to an adiabatically constant value. Also, as ex- 
pected, the type III mode exhibits an oscillating behavior 
before entering the resonance band. 

We can also conclude from Figure 5 that the particle 



number follows the predicted evolution inside the reso- 
nance band, the main differences between the numerical 
and the analytical solutions occurring near the endpoints 
of the resonance band, where the transitions between res- 
onant and oscillating regimes take place. The discrep- 
ancies for low particle number are also due to the nor- 
malization of rife, as discussed earlier. It is nevertheless 
clear that Eq. (|100p gives a quite good description of 
the resonant particle production regime in an expanding 
universe, although it underestimates the particle number 
density, and that the contribution of the non-resonant 
periods are subdominant, as expected. 

The total particle number produced by the end of the 
resonant regime, (22), was computed for several values 
of /3, keeping |qo| and 7 fixed at the values chosen above. 
The results we have obtained are shown in Figure 6. 
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W The positive value of log reflects the normalization chosen in 
Eq. II103H and does not correspond to a physical particle number. 
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FIG. 6: Numerical results for the value of the total particle 
number produced by the resonance as a function of f3. Plot 
(a) shows the results for small values of (5, including type 
II and III modes, and plot (b) shows the results for large 
values of /3, where only type III modes are present. The solid 
line corresponds to the numerical results while the dashed 
line gives the corresponding analytical prediction. All results 
correspond to |qo| = 10~ 3 and 7 = 1CP 7 . 

One observes that hk (^2) increases with [3 for type II 
modes and slowly decreases with (3 for type III modes, 
in agreement with Eq. (|101[) . The main discrepancies 
between the numerical results and the analytical predic- 
tion are again due to the smooth transition between the 
resonant and non-resonant regimes at Z2, with the latter 
giving a subdominant contribution to the total particle 
number, which is more significant for type III modes. The 
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different number of oscillations the modes undergo before 
entering the resonance band give the oscillations that can 
be observed in Figure 6. These are suppressed for large 
[3 as the pre-resonance oscillations of the particle number 
give a negligible contribution for high momentum modes. 

Numerical simulations also show that, as expected, the 
particle number in all modes decreases with 7, so that 
Hubble expansion may inhibit the resonant amplification 
of particle modes at early times. 

Thus, we conclude that our analysis of the resonance 
regime gives a good description of the mechanism of par- 
ticle production in branonium for an expanding universe 
and that the corrections arising from the non-resonant 
regimes are subdominant. Extrapolation of this discus- 
sion to high momentum modes is, however, difficult, as 
these will only be excited at late times, where the approx- 
imations we have considered may no longer hold. The 
production of high momentum particles is nevertheless 
suppressed, as we have concluded earlier. Furthermore, 
at each given time z > 0, only the modes with z\ < z 
have entered the resonance band, which means that only 
those modes for which 



Aok < A Qk 



i + |«b|(i + £s) 
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have started being excited. This gives the following mo- 
mentum cut-off for modes which have already been am- 
plified at t > t : 
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(105) 

Notice that this gives the expected endpoint of the reso- 
nance band in the limit a — > 0, given by Eq. (|73|) . 

The generalization of these results for massive parti- 
cles follows the same reasoning used in the non-expanding 



case, with the substitution k 2 hys 



kp hys + m 2 within our 
the 



approximations. Note that, in the expanding case 
physical momentum of the particles is redshifted while 
their mass remains the same as the universe expands. 
This modification will then produce a more complex evo- 
lution of the parameter Ak(z) in the modified Mathicu 
equation, Eq. (|9ip . given to lowest order by: 



string values. We need, however, to bear in mind that 
Eq. (|91[) is not the Mathieu equation and that several 
other terms in Eq. (|52| may become important at early 
times. These may alter the excitation time and the am- 
plification of each mode, or even prevent any particle 
production. It is nevertheless clear that the effects of 
the expansion are suppressed at the typical energies in- 
volved in the probe's motion, so that one still expects a 
significant number of particles to be produced. 



C. Energy radiated into brane particle- modes 

As the resonance develops, the probe brane loses en- 
ergy to excited 77-particle modes and one may wonder 
whether the energy radiated by this mechanism is suffi- 
cient to affect the probe brane's motion. 

Let us start by computing the energy density of 
massless particles produced by the resonance in a non- 
expanding universe. We have seen that each rescaled 
mode function Xk has an associated harmonic oscilla- 
tor energy function E k — ^fe(^fe + 1/2). For the mode 
functions Xfc = \ffX k , the associated energy function is 
approximately the same in the non-relativistic and large 
distance limit, where / ~ 1. Hence, the total energy 
density radiated into brane particle-modes is given by: 



d 3 k 1 
J2wf k ~ 2^2 
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dkk z cj k [n k + - ) , (107) 



where we have used the momentum space isotropy dis- 
cussed earlier in this work. After a significant number 
of particles has been produced, we may take n k 3> 1/2 
and neglect the 1/2 factor in Eq. (I107j) . Also, we know 
that particles are produced in a narrow resonance band 
centered at k c = fi/2 with a band width Ak res ~ <5 2 /fi = 
|g|fi/2. The particle number distribution in momentum 
space can then be well approximated by a Gaussian dis- 
tribution of the form: 



. 1 (fc fcc) 
nfc (i)=n fce (t)exp( -2^7^ 



(108) 
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where s = ^7 and all other quantities are defined as 
before. Although this will alter the way the resonance 
band moves into higher momentum modes as the probe's 
orbit decays, the evolution of the modes is expected to 
follow the same qualitative behavior as in the massless 
case, so that we will not analyze this in more detail. 

We, thus, conclude that, to lowest order, the paramet- 
ric resonance survives in an expanding universe, although 
modes are excited only during a finite period. These re- 
sults hold only at late times but we note once again that 
all quantities are measured with respect to the typical 



The gaussian distribution width Afc should correspond 
to a fraction of the resonance band width, so that we 
may write Afc = XAk res . The coefficient A is expected to 
be time-dependent as the central mode fc c is more ampli- 
fied than all the other modes in the resonance band, so 
that the distribution should become more sharply peaked 
about this value with time. The advantage of writing the 
particle number distribution in this form is that it allows 
us to further approximate the result by a <5-function dis- 
tribution, taking into account that Afc/fc c <C 1, as dis- 
cussed earlier. Introducing the correct normalization, we 
may then write: 



(109) 



n k {t) ~ n kc (t)v / 2TrAk(t)S(k - k c ) . 
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Using coke — k c — 0/2, we obtain: 



(110) 



where the time dependence of nfe c , A and consequently 
e p is implicit. Notice that, when q — > 0, we have e p — > 
0, as the resonance vanishes in this limit. To have an 
estimate of the order of magnitude of this quantity, take 
the example we have considered before in the discussion 
of particle production in a non-expanding universe, for 
which \q\ — 10~ 4 and fl — 4 x 10~ 5 , in units of the string 
length. By solving the equations of motion numerically 
and fitting the obtained particle number distribution to 
the gaussian distribution given in Eq. (|108[) . we may 
determine the values of A and nk c . For an initial time 
to = 10 9 , we obtained the following results: 



n kc (Wt ) 



7.5 x 10 2 

1.6 x 10 7 



A(5t ) : 
A(10io) 



0.20 , 

- 0.12 . (Ill) 



These two examples confirm that, indeed, the reso- 
nance width decreases in time. From Eq. (|110|) . we 
obtain for these two cases e p (5to) ~ 3 x 10 -22 and 
e(10fo) — 4 x 10~ 18 , in units of the string energy scale. 

The energy density radiated into brane particle-modes 
should then be compared with the energy density of 
the interbrane distance field cj>, which is given by e = 
—a/(2R) = — (crSl/2) a . For the example considered 
above, e ~ —2 x 10 -2 , so that the value obtained for the 
energy density radiated into particles up to t — 10£o = 
10 10 is negligible when compared to the probe brane's 
energy density. In general, we obtain: 
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which gives the main result of this subsection. This 
means that, in the previous example, if we take A ~ 
0.01 — 0.1, up to a ^-particle number density of 10 23 — 10 24 
can be produced without affecting significantly the mo- 
tion of the probe brane. 

These results can be easily generalized for the massive 
case, by taking into account the changes in k c and Ak res 
discussed earlier. For S 2 /2 <C fc 2 , we obtain: 
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with k r 



and which, as expected, re- 



duces to Eq. (|112|) in the limit m — > 0. As the momen- 
tum of the produced particles decreases with m, it is easy 
to conclude that the energy damped into brane particle- 
modes is smaller for more massive particles (recall that 
the strength of the resonance is not affected by the mass 
of the particles). 



The analysis of the energy radiated into brane particle- 
modes is more difficult to perform in the expanding uni- 
verse case, as more effects have to be taken into account. 
First, we need to recall the approximate relation between 
hk and the physical particle number. Next, we need to 
notice that the mode functions Xk are redshifted by a 
factor of <z(i)~2 with respect to the harmonic oscillator 
mode functions Xk- Finally, for the purposes of deter- 
mining the energy density of the produced particles, we 
may take u>k(t) ~ k/a(t). Then, we have 
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We see that the energy density of the produced parti- 
cles is redshifted by the usual power of a -3 , but is further 
reduced due to the redshift of the physical momentum 
and energy of the modes. We can use the expressions ob- 
tained for the particle number produced by the resonance 
in Eq. (|10ip to determine hk, as we have concluded that 
the effects of the non- resonant periods are subdominant. 
Writing them in terms of the comoving momentum of the 
modes, we conclude that, for type III modes, 
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To determine the limits of this integral, one must re- 
call that, at a given time, only a finite number of modes 
has entered the resonance band. In particular, one may 
use the value computed in Eq. (| 1 05[) as the upper limit 
of this integral, although some of these modes have not 
yet been completely excited by the resonance. The lower 
limit of the integral will be the first type III mode, with 
k = fc* = fccoyl + l<7o|- This gives a finite energy den- 
sity radiated into brane modes at each finite time t > tg. 
This energy density is redshifted by Hubble expansion, 
but at the same time modes with increasing energy en- 
ter the resonance band. This process will, however, stop 
at some point, as our approximations will break down 
when the motion of the probe becomes relativistic and 
gets too close to the central stack. Even if the resonance 
mechanism persists in this limit, although certainly in a 
different form, the probe will eventually modify its tra- 
jectory, either by colliding with the central stack or via 
some other mechanism that stabilizes its motion. 

The contribution of type II modes is obviously finite 
and can be computed by integrating Eq. (| 1 14[) using the 
corresponding expression for hk given in Eq. (jlOip , with 
a lower limit k c o and an upper limit fc* for the integral. 
No particular insight is, however, gained by computing 
the exact expressions for the energy density for both type 
II and type III modes. Instead, we note that the domi- 
nant contribution to the energy density will be given by 
high momentum modes which enter the resonance band 
at late times. For these, the effects of the expansion of 
the universe will be less significant and one may use the 
results obtained in the non-expanding case to compute 
the energy density produced by the resonance. In this 
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case one should use the values of q and fl at the end of 
the resonance mechanism. 



D. Radiation into bulk modes 

As discussed in [l(J, the probe's trajectory in the back- 
ground created by the central stack may induce not only 
the production of brane-bound particles but also radi- 
ation into bulk modes, namely gravitational, RR-form 
and dilatonic fields. This effect is due to the acceler- 
ated motion of the probe brane and most of the power is 
radiated into lower-spin fields, in particular scalar fields 
which couple to the brane's orbital monopole moment. 

The probe will thus lose energy through this process, 
contributing to the decay of its orbit. This will certainly 
modify the mechanism of resonant production of particle- 
modes in the brane, namely by varying the orbital fre- 
quency, which determines the comoving momentum of 
the modes that are excited by the parametric resonance. 
One then expects the effects of radiation into bulk modes 
to be quite similar to those induced by the universe's ex- 
pansion, so that we will not analyze this process in detail. 
Nevertheless, it is important to estimate its contribution 
to the decay of the probe's orbit. 

The power radiated into bulk scalar modes was esti- 
mated in [10], where it was shown that the number of 
orbits the probe can complete before the interbrane dis- 
tance becomes of order l s is approximately given by: 

at- AteV—L-, (He) 

where is the initial value of the physical radius of the 
orbit, assuming it is circular. (This gives a good esti- 
mate for the decay time, even though we are interested 
in small eccentricity orbits for the resonance process. Af 
corresponds to the ratio of the decay time to the ini- 
tial orbital period.) This expression shows that, for the 
probe to complete a large number of orbits before decay, 
it needs to be at a large distance from the central stack, 
in units of the string length. Also, one needs c/ 3 N <C 1 
to obtain a sufficiently large value of Af. 

This value should be compared to the number of orbits 
of the probe necessary for the resonance to be effective, 
which in the non-expanding case is of order (2ir^,k)^ 1 ~ 
(tt I Qf I ) 1 . Hence, the production of a significant number 
of particles requires Af ^> (7r|g|) _1 , which gives for the 
initial radius of the orbit: 

Ma)'"'- (117) 

For example, if |g| ~ 10~ 3 and there are 10 branes in 
the central stack, we need the initial interbrane distance 
to be larger than about 300<? s ^ s for the resonance to be 
effective, which is not too large a number taking into 
account that g s is parametrically small. 



As, in the analysis of the parametric resonance mech- 
anism, we have assumed the interbrane distance to be 
large compared to the string length, we expect the ef- 
fects of radiation into bulk modes to be initially negligi- 
ble. This process should, however, become more impor- 
tant as the orbit decays and increases its acceleration. It 
is even possible that, at late times, it overcomes Hubble 
expansion as the main energy loss process. Nevertheless, 
the end of the resonance should still be determined by 
the break down of the non-relativistic and large distance 
approximation, as this also controls the amount of energy 
damped into bulk modes. 

IV. EFFECTS OF TRANSVERSE SPACE 
COMPACTIFICATION 

So far we have considered the probe to move at dis- 
tances from the central stack which are small compared 
to the typical size of the transverse space directions, so 
that we may neglect the effects of compactification on 
the probe's motion. In this section, we will analyze the 
leading order effects introduced by the finite size of these 
directions, and show that they may lead to the creation 
of orbital angular momentum, a necessary condition for 
resonant particle production to take place. 

We will consider the simplest case of compactifying the 
three transverse directions on an isotropic 3-torus of size 
R± , which is defined by the identifications y l <-* y l + R± . 
As mentioned earlier, the harmonic function associated 
with the central stack configuration needs to be modi- 
fied in this case by including the appropriate "brane im- 
ages", according to Eq. (fTU|) . These images are placed 
at points with coordinates y l = riiR±, with integer m, 
defining a hypercubic lattice corresponding to the cov- 
ering space of the 3-torus, (R/Z) 3 . Although toroidal 
compactifications are too simple to give realistic particle 
physics phenomenology, this will be sufficient to illustrate 
the main effects of compactification on the probe brane's 
trajectory. 

The interbrane potential associated with the general- 
ized harmonic function Eq. (fTt?|) is obtained by comput- 
ing the propagator for a massless field on a torus, which 
satisfies: 

V 2 G(y,y')^(y-y')-T^ , (118) 

where V± = R^ is the volume of the transverse 3-torus. 
The term —1/V± is included for consistency, so that 
the integral over the compact manifold of both sides of 
Eq. (|118p vanishes. If one expands the Green's function 
G(y, y') in terms of eigenfunctions of the Laplacian oper- 
ator V 2 , one concludes that the — l/V± term removes the 
unphysical zero-mode which makes the Green's function 
diverge. This term arises naturally from the curvature 
of the non-compact 4-dimensional spacetime [24| . It has 
been shown that one can write the massless propaga- 
tor on the 3-torus in terms of an integral involving the 
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Jacobi-theta function 9 3 



G(y) = G(y 



(119) 

Although this expression gives the full compact space 
propagator, we are mainly interested in the leading or- 
der modifications introduced by compactification when 
the probe is moving closer to the central stack than to 
any of its images (orbits around any of the brane images 
are equivalent due to the symmetries of the hypercubic 
lattice). In (25[, an expansion of this propagator about 
the origin was computed using Ewald's method for calcu- 
lating potentials in hypercubic lattices in the context of 
solid state physics. The leading corrections to the central 
1 jr potential are given by: 
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The values for the coefficients of the various terms were 
found numerically to be C$ = —0.21, A 4 = 0.44 and 
A e = 0.0072. The functions h 4 (y) and h e (y) are har- 
monic functions of order (r/R±) 4 and (r/R±) 6 , respec- 
tively, and the series continues with harmonic terms of 
higher (even) orders. The leading harmonic correction to 
the propagator is given by: 
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Recalling that the function h(y) contributes to the 
interbrane potential through the graviton-dilaton and 
RR-form interactions, we may write this potential for 
Qe <C r <C R_l, to order (r/R±) 4 and discarding con- 



stant terms, as 
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where we defined the function 
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(123) 



We conclude that, to this order, the interbrane poten- 
tial is modified by two terms. Both terms correspond to 
repulsive contributions to the potential arising from the 
overall attraction of the image branes and their coeffi- 
cients are suppressed by a factor 1 / R± , so that their ef- 
fects are negligible for large transverse volume. The "jel- 
lium" term, of order (r/R±) 2 and whose name arises in 
the context of solid state physics, gives an isotropic con- 
tribution while the "asymmetry" term, of order (r/R±) 4 , 
gives an anisotropic contribution to the potential. 



We will be more interested in the effects of the asym- 
metry term, as its breaks the rotational symmetry of the 
transverse space and leads to the generation of orbital 
angular momentum. The effects of the jellium term will, 
however, be more significant close to the central stack. 
It is easy to show that, for large transverse volume, this 
term will, to leading order, make the probe brane's or- 
bits precess. This is due to its repulsive nature and, 
for small eccentricity orbits, the deficit angle is approx- 
imately given by A9 ~ 32tt 2 (Z 2 / R±Q 6 (T 6 V 3 ) 2 ) 3 . This 
effect will, however, be suppressed at late times as the 
probe's angular momentum is rcdshiftcd away by Hub- 
ble expansion. 

The inclusion of the asymmetry term makes the anal- 
ysis quite difficult, as the potential is no longer central 
and depends on both angular coordinates in the trans- 
verse space. This may lead, in particular, to non-planar 
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orbits of the probe brane. However, if we set y 
initially, the probe will feel no force along this direction 
and its trajectory will be confined to the (y , y 2 ) plane as 
before. We will focus on this particular case, bearing in 
mind that in general non-planar trajectories may arise. 

Within these assumptions, the problem reduces, as be- 
fore, to the evolution of the canonically normalized com- 
plex scalar field 4>, defined in Eq. (jTTJ) , in an expanding 
universe. Its potential can be written as 
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where cf>R and <pi are, respectively, the real and imaginary 
parts of the field, a was defined in Eq. (|2T|) and 
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give the strength of the jellium and asymmetry terms. 
This potential is no longer invariant under the global 
U(l) symmetry of the non-compact case, signaling that 
the probe's angular momentum is no longer conserved. 
This is quite similar to the U(l)-violating potential gov- 
erning the evolution of the scalar field in the Affleck- 
Dine mechanism (26j , where non-conservation of angular 
momentum plays a crucial role in generating the U(l) 
baryon number asymmetry in our universe. The po- 
tential remains, however, invariant under the hypercubic 
group symmetries 4>rj — > ~4>r.i and 4>r <-></>/. In terms 
of the polar angle 9, this means that it/ 2 rotations as 
well as reflections about the 9 = ir/4 axis are preserved 
by the asymmetry term, as one can easily conclude by 
writing Eq. (|123[) for y 3 = as 



f(0) = 1 - 2 sin 2 (2(9) 



(126) 



Hence, we only need to consider initial conditions such 
that <f)R > <f)i > or equivalently < 9 < n/4. The 
U(l) symmetry violation precludes a complete analytical 
description of the orbits. However, it is not difficult to 
obtain a qualitative insight on the main features of the 
probe's trajectories. If the probe is placed significantly 
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far from the central stack (but still closer to it than to 
any of the image branes), the asymmetry term will not be 
negligible and some angular momentum will be created 
or destroyed. If the probe has no angular momentum 
initially, as we would expect immediately after inflation, 
then it will necessarily acquire some and be placed in an 
orbit around the central stack, instead of just falling to- 
wards it along the radial direction. As both the jellium 
and the asymmetry interactions are repulsive, it is possi- 
ble that this trajectory does not remain bound to the cen- 
tral stack, becoming connected to one or more branes in 
the hypercubic lattice. Hubble expansion will, however, 
redshift the probe's energy and angular momentum, so 
that one expects the probe to become bound to only one 
brane stack at late times. As the orbital radius decreases, 
the U(l)-violating term becomes less significant and the 
particle number associated with <f> should asymptotically 
become constant, i.e. the angular momentum should vary 
only due to the universe's expansion. One also expects 
the jellium term, as well as the asymmetry term, to in- 
duce some precession of the orbital axis, which should 
stop at late times when these terms become negligible. 

Due to the hypercubic symmetries, along 6 = 0,7r/4 
the force acting on the probe will be in the radial direc- 
tion, as /'(7r/4) = /'(0) = 0, and no angular momen- 
tum will be created in these particular directions. For 
< 6 < 7r/4, angular momentum creation should, how- 
ever, be a generic feature. 

To have a better understanding of how the initial con- 
ditions affect the amount of angular momentum created 
or destroyed, we have computed the force acting on the 
probe brane and the associated torque, given by: 
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where we have considered the field as a vector in the 
(<Pr, 4>i) plane. Recall that the torque gives the variation 
of the angular momentum, which in terms of the particle 
number density n defined in Eq. (|20[) can be written as: 
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We have plotted the isotorque contours for some par- 
ticular parameter values in Figure 7. We have restricted 
the position of the probe to the region y 1 ,y 2 < 0.5i?^, 
where the expansion of the potential about r = is suffi- 
ciently accurate and the probe's motion is bound to the 
central stack (near the limits of this region, the orbits 
may not be bound to the central stack, as the actual re- 
gion where the force points towards the central stack has 
a more complicated shape). 

Observing this figure we see that the torque is max- 
imized close to the boundary of the region mentioned 
above (darker areas for < 8 < ir/4). Also, the longer 
the probe's trajectory remains close to this boundary the 






FIG. 7: The torque acting on the probe brane in the (y 1 ,y 2 ) 
plane. Coordinates are normalized to the size of the trans- 
verse torus, with the value R± = 2 x 10 4 in this example, 
in units of the string length. We considered Q§ = 100 and 
TeV3 = 1 in these units. Note that the torque changes sign 
when crossing one of the 9 — nn/4 axis, with integer n, being 
negative for < 8 < 7r/4. Also, the absolute value of the 
torque increases as one moves away from the origin. 



more angular momentum it is likely to gain. Thus, max- 
imum angular momentum creation should occur if the 
probe is placed below the centre of the hypercubic cell. 

The initial value of the Hubble parameter H will also 
affect the amount of angular momentum created. If it is 
of the same order of magnitude as the orbital frequency, 
one expects the probe to be driven into lower torque re- 
gions within a few periods. Lower values of the Hubble 
parameter should allow the probe's angular momentum 
to oscillate significantly as the probe moves through al- 
ternate regions of positive and negative torque, but no 
net angular momentum will be gained after completing 
the first few orbits. The asymptotic value of the comov- 
ing particle number N — na 3 will, hence, be determined 
by the initial path of the probe. 

We have simulated the evolution of the probe brane nu- 
merically to illustrate this discussion. The initial values 
for the real and imaginary parts of the field were defined 
as follows: 
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with 8r > Si, giving < 9 < 7r/4, and Srj < 0.5, which 
ensures the probe is placed initially within the region 
where it is bound to the central stack. We have also 
computed the eccentricity of the orbit, as defined in Eq. 
([29]) . with the probe's energy density including both the 
jellium and asymmetry terms. This quantity should only 
become meaningful at late times, when the latter terms 
become negligible, but its evolution tracks the creation 
of angular momentum along the probe's trajectory. 

Figure 8 shows an example of the results obtained for 
the probe's orbit. In this case, the probe is placed be- 
low the centre of the hypercubic cell and, as expected, 
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FIG. 8: Results obtained numerically for the probe brane's 
trajectory for Q 6 = 100, R± = 2 x 10 4 and T 6 = Vf 1 = (2tt) 6 . 
The probe is initially stationary at the point (Sr,5i) — 
(0.50,0.35) for to = 2 x 10 6 and we consider a matter- 
dominated universe with a = 2/3. The asymptotic value 
obtained for the eccentricity is e ~ 0.28. 



a significant angular momentum is produced during its 
motion. One observes, as expected, a small precession of 
the orbital axis and the decay of the orbit into regions 
of lower torque, which drives the orbital eccentricity to- 
wards a constant value. Initial conditions were chosen in 
this example so that the inverse Hubble parameter is of 
the same order of magnitude as the orbital period and, 
hence, the orbit is quickly stabilized. The eccentricity 
tends in this example to a constant value of e ~ 0.28. 
Notice the similarities with the example illustrated in 
Figure 1, the main difference residing in the fact that, in 
this case, the probe has initially no angular momentum. 

In Figure 9 we plot the asymptotic values of the eccen- 
tricity obtained numerically for different initial positions 
in the plane (<f)R,<f>i), considering two distinct values of 
the initial Hubble parameter in a matter-dominated uni- 
verse. In both cases, one observes that smaller eccentric- 
ity orbits are obtained near the larger torque regions be- 
low the centre of the hypercubic cell, as expected. We see 
that, in general, some angular momentum is produced, 
although the majority of initial conditions leads to highly 
eccentric orbits, with e > 0.9. We also conclude that 
the probe's orbit is generically less eccentric if its motion 
starts at later times, so that the initial Hubble damping 
of the asymmetry term is smaller. As one may observe in 
Figure 9, in particular in plot (b), the combination of the 
asymmetry term and Hubble expansion makes the final 
eccentricity vary in a non-trivial way with the initial posi- 
tion of the probe, namely near Sr ~ 0.5, where the orbits 
may in some cases become unbound. Note that these re- 
sults correspond to particular values of the parameters of 
the problem and that the values of the asymptotic eccen- 
tricity are highly dependent on these parameters, namely 
the relevant distance scales Qq and R±. 
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FIG. 9: Contour plots obtained numerically for the asymp- 
totic value of the orbital eccentricity in the plane (Sr, Sr), 
with Sr > 8 1 (with similar results for Sr < Si), assuming 
Q 6 = 100, R± = 2 x 10 4 and T 6 = l/f 1 = (2tt) 6 . Both plots 
were obtained for a matter-dominated universe with a = 2/3 
and the motion of the probe starts at (a) to — 1 x 10 6 and (b) 
to — 5 x 10 6 . All parameters are given in units of the string 
length. Contours were plotted for e = n/10, with n = 3, . . . ,9, 
and darker regions correspond to lower eccentricities. 



These results show that a significant amount of angular 
momentum can be created by U(l)-violating terms in 
the probe's potential arising from compactification of the 
transverse dimensions, which is a key ingredient for the 
development of the parametric resonance we have studied 
in this work. The simple 3-torus compactification thus 
ilustrates how low eccentricity orbits can be created for a 
range of initial conditions, even if the probe has initially 
no angular momentum. 



V. STABILIZATION AND COSMOLOGICAL 
IMPLICATIONS 

Potential applications of the resonant particle produc- 
tion mechanism to cosmology will be largely determined 
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by the system's final configuration. As we have discussed 
before, the D6 — D6 branonium we have analyzed in this 
work is unstable, as the probe (anti)brane loses energy 
and angular momentum due to Hubble expansion. At 
some stage, the distance between the probe and the cen- 
tral stack of branes will consequently become of the or- 
der of the string length (although our approximations will 
break down long before this happens) . When this occurs, 
the scalar mode associated with open strings stretching 
between the probe and the source branes becomes tachy- 
onic, signaling an instability of the system (28[. Con- 
densation of this tachyonic mode will then lead to the 
annihilation of the probe brane with one of the branes 
in the central stack. The non-supersymmetric nature of 
the D6 — D6 configuration, which is a general property of 
all systems mixing brane and antibrane states, is behind 
this instability. Namely, it makes the system evolve into 
a 1/2 BPS state with (JV — 1) parallel DQ branes, with 
the same total charge as the original configuration but 
with a lower energy. 

Interesting scenarios may arise if some additional 
mechanism stabilizes the probe brane at a finite distance 
from the central stack, in which case the produced parti- 
cles may survive after the resonance mechanism ends. An 
attractive possibility for stabilizing multiple brane sys- 
tems was suggested in [27j | . When supersymmetry break- 
ing occurs, one or more fields within the supermultiplets 
of the theory typically acquire masses, their value being 
set by the energy scale at which SUSY is broken. Massive 
gauge potentials will then induce short-range interactions 
which, if repulsive, may balance the gravitational atrac- 
tion between different branes. Notice that the presence 
of a probe antibrane itself breaks SUSY but that this 
does not induce bulk field masses, so that other sources 
of SUSY breaking need to be considered in this case. 

Suppose, for example, that the C\ RR-form and conse- 
quently its magnetic dual form CV become massive after 
SUSY breaking, the same happening to the dilaton field, 
$. Both the RR-form and dilaton- mediated interactions 
will then be described by Yukawa potentials whose range 
is determined by the masses mun and m$, respectively. 
If turr <C m$, the dilaton- mediated interaction will be 
exponentially suppressed for interbrane distances of or- 
der m~^ R . Assuming Qq <§; rn^ R <C R±, we may write 
the relevant terms in the potential for r • 
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as 



V{r) = M A 
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where M 4 gives the overall constant factor. The con- 
stant A parametrizes the effective strength of the RR- 
interaction compared to the gravitational atraction and 
one expects |A| > 1 as the dilaton- mediated interaction 
is negligible. If, as we have assumed so far, the probe is a 
£>6-brane, the RR-interaction, although short-ranged, is 
attractive and cannot balance the 1/r gravitational part. 
If, however, the probe is a _D6-brane, we have A > and 
stabilization of the system may be possible. In this case, 



for r ^> m RR , the RR-interaction is exponentially sup- 
pressed and the probe's motion is governed only by the 
gravitational term. The evolution of the probe brane will, 
in this case, be very similar to that of a probe antibrane in 
a massless RR-form background which we have studied in 
this work, the potential having a smaller strength than 
in the latter case. This difference does not modify the 
qualitative features of the probe's trajectory and all the 
results we have derived in this work remain valid using an 
effective central stack charge Q e J^ < Qe- In particular, 
one expects the parametric resonance to develop in the 
probe's world-volume under the same conditions as be- 
fore and the effects of transverse space compactification 
to provide the necessary source of angular momentum for 
similar initial conditions. 

The fate of a probe Z?6-brane will, however, be quite 
different than that of a probe £)6-brane. As the probe's 
orbit decays due to the universe's expansion, the inter- 
brane distance will eventually become of order rrfj^. At 
these distances, the repulsive RR-interaction becomes 
relevant to the probe's motion and stabilization may be 
possible. The potential given in Eq. (|130[) has a local 
minimum formally given by: 
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where W(z) is the Lambert W-function, defined as the in- 
verse of f(W) = We w . After its angular momentum has 
been redshifted away by Hubble expansion, one expects 
the probe to settle at this local minimum. This may not 
be the absolute minimum of the full potential, which in- 
cludes the dilaton-mediated interaction and other super- 
gravity and string theory corrections closer to the central 
stack. Hence, absolute stabilization of the probe cannot 
be guaranteed, but it is reasonable to expect at least a 
long-lived metastable state. 

From the properties of the Lambert W-function, one 
concludes that ro = for A = 1, being positive for A > 1, 
and that ro strictly grows with A. Hence, as expected, 
stabilization at large distances from the central stack is 
only possible if the RR-repulsion is stronger than the 
gravitational attraction. In this case, one expects the 
probe brane to stabilize away from the central stack at a 
distance ro ~ m RR- 

The scale of supersymmetry breaking will then deter- 
mine the cosmological implications of the particle pro- 
duction mechanism analyzed in this work. If SUSY is 
broken at very high energies, inducing a large mass for 
the RR-potential (but still smaller than the dilaton's 
mass), the probe will stabilize very close to the central 
stack and no tachyonic modes are present to induce an- 
nihilation. The stable probe and the central stack will 
then be characterized by a general broken gauge symme- 
try group U(N) x U(M), where M < JV is the number 
of branes in the probe, with an exponentially large num- 
ber of particles charged under the U(M) gauge group if 
the resonance develops before stabilization occurs. This 
may be relevant, for example, to the generation of the 
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baryon asymmetry in our universe. On the other hand, 
a parametrically small soft SUSY breaking mass will sta- 
bilize the probe at large distances. The parametric reso- 
nance may then be relevant for dark matter production, 
as particles in the probe's world-volume will necessarily 
interact weakly with the visible sector if the latter is em- 
bedded in the central stack. As we have shown, these 
particles may be supermassive (although parametrically 
small compared to the string scale) , so that this scenario 
could provide an effective mechanism for a non-thermal 
production of WIMPZILLAS H- 

One also expects the parametric resonance to amplify 
other bosonic modes living in the probe's world-volume, 
such as Yang-Mills fields if the probe has more than one 
brane. If these mediate baryon number-violating forces, 
their interactions with central stack fields could then pro- 
vide a mechanism for baryogenesis as the one discussed in 
[3(il ] , in which case the probe would need to be stabilized 
before nucleosynthesis. 

It is also possible for fermionic particles living in the 
probe's world-volume to be produced in resonance, al- 
though Pauli blocking makes this process significantly 
different from the bosonic case, according to the discus- 
sion given in [3l|. Thus, if the probe can be stabilized, 
a resonant production of fermionic dark matter particles 
may also be achieved through this process. 

One also expects the interbrane distance field to oscil- 
late about the local minimum of the potential, possibly 
generating a second stage of resonant amplification. This 
may be induced from coupling to other fields, as in the 
standard preheating mechanism, or gravitationally, as in 
the case described in this work. Although this second 
stage may contribute significantly to the final particle 
number density in the probe brane's world-volume, we 
will not discuss them in further detail but rather refer 
the reader to the discussions given in this work and in 
the literature [Hd03. 

We emphasize that all these scenarios depend on the 
particular embedding of the Standard Model fields in this 
setup, an issue for which there is still no complete answer 
and that is closely related to the nature of SUSY breaking 
and its communication to the Standard Model. 



VI. CONCLUSIONS 

In this work we have analyzed the mechanism of par- 
ticle production in D6 — D6 branonium systems, as well 
as in D6 — D6 systems with SUSY breaking. We have 
shown that, in the limit of large distances and small ve- 
locities and for small eccentricity orbits, a parametric 
resonance will develop, producing scalar particles con- 
fined to the probe brane's world-volume. Massless par- 
ticles are produced in a narrow resonance band centered 
at the comoving momentum corresponding to one half 
of the probe's orbital frequency f2, as if resulting from 
the direct decay of the particles associated with the in- 
terbrane distance field. The strength of the resonance is 



given by the exponent /ifc defined in Eq. (|72|) . Massive 
particles are produced in the same way, with energies of 
order tt/2. These energies are small compared to the 
string energy scale but, as the latter should be close to 
the Planck scale, very massive particles can be produced 
by the resonance mechanism, according to Eq. (|84| . For 
both massless and massive particles, the associated parti- 
cle number grows exponentially in time and a large num- 
ber of particles can be produced by this mechanism, just 
as in the case of preheating after inflation [H, El ■ How- 
ever, the resonance surely requires the probe to complete 
a large number of orbits before the particle-modes can 
be significantly amplified, due to the small value of the 
exponent /x^ ~ § <C 1. 

Hubble expansion of the 4-dimensional effective flat 
FRW spacetime makes the probe's orbit decay, redshift- 
ing its energy and angular momentum. This alters the 
development of the parametric resonance significantly if 
the probe's motion does not occur at sufficiently late 
times. If, however, its motion begins late in the his- 
tory of the universe, at least compared to typical string 
times, the resonance will still develop with the associated 
resonance band being shifted towards higher momentum 
modes. Each mode will be excited during the finite pe- 
riod it spends inside the resonance band, after which the 
associated particle number becomes an adiabatic invari- 
ant. A significant number of particles, given appproxi- 
mately by Eq. (|I0I[) . can be produced if this period is 
sufficiently long for the resonance to develop, according 
to the condition > 7, where 7 quantifies the effects of 
the expansion and was defined in Eq. (|93[) . 

We have also concluded that a large number of brane- 
bound particles can be produced without affecting the 
probe's orbital motion, as can be seen in Eq. (|II2p . Ra- 
diation into bulk closed string modes will also have a 
negligible effect if the probe moves at sufficiently large 
distances, as assumed for the validity of our study. 

A realistic implementation of this mechanism requires 
particular initial configurations with large interbrane dis- 
tances, non-relativistic velocities and almost circular or- 
bits. We have shown that angular momentum creation 
may result from the effects of compactification of the di- 
rections transverse to the branes, illustrating this for the 
case of a compact 3-torus, where low eccentricity orbits 
can be produced for some range of initial conditions. 

In this work, we have also discussed how, for a probe 
_D6-brane, an interbrane potential may arise from super- 
symmetry breaking, if e.g. the relevant RR-form and 
the dilaton, <!>, gain a mass in the process. Stabilization 
of the probe is then possible at late times and at dis- 
tances of order m^ R , if m$ ^> mun- Such a stabilization 
avoids annihilation with the central stack, as necessarily 
happens for a probe antibrane, making the parametric 
resonance regime potentially relevant for baryon number 
or dark matter generation, for example. 

In our analysis, the parametric resonance arises only 
for D6 — D6 or D6 — D6 systems, where the 1/r poten- 
tial allows the probe to move in closed elliptical orbits in 
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a non-expanding universe. It is crucial for the develop- 
ment of the resonance that the interbrane distance field 
p exhibits an oscillating behavior, although the orbits do 
not need to close, as we have seen for the case of an 
expanding universe. Thus, other Dp — Dp or Dp — Dp 
branonium systems with p < 6, where the probe trajec- 
tories are necessarily unbounded from below, cannot ex- 
hibit such resonant particle production. Additional sta- 
bilization potentials may, however, provide the required 
oscillatory behavior, so that it may be possible to find a 
resonant particle production mechanism in other cases. 

It is nevertheless clear that resonant particle produc- 
tion in branonium systems may play an important role 
in our universe's evolution and we hope with this work 
to motivate future exploitation of their properties and 
applications. 
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APPENDIX A: COMPUTATION OF THE 
PARTICLE NUMBER FOR H >0 

The particle number produced by the parametric res- 
onance can be estimated by computing the integral in 
Eq. (|100p . This gives the physical particle number apart 
from constant factors, subdominant non-resonant contri- 
butions and also an a~ 3 redshift factor due to Hubble 
expansion. We have, 



{0, z < Zi , 

2 J^fj, k (z')dz', Zl <z<z 2 , (Al) 
2 J z . 2 fi k (z')dz', z > z 2 ■ 

In order to compute this function, let us write the ex- 
ponent /i k (z) in the form: 



where 



a 
b 

c 



<&e - aw , 

2(<?o 2 £ + A kQ (A k0 - 1) 7 ) , 
q 2 - (A k0 - If . 



(A3) 



For all modes of interest, a < 0, while we have c > for 
type II modes and c < for type III modes. The sign 
of b will not affect our results. It is also useful to define 
d 2 = b 2 — Aac, which is explicitly given by 

d 2 = 4q 2 (2A k0 (A k0 -l)^+(A k0 -l) 2 e+A 2 0l 2 ) , (A4) 

and is positive for all type II and type III momentum 
modes. With these considerations, we find, for Z{ < z < 

Z2, 



lognjTW 



2<i 



(2az' + b)/j, k (z') + 



<f : arcsin f 2 ^ 3 ±I> 



(A5) 



Taking into account that fj- k (0) — \fc/2, p k (zi) = 
Hk{z2) = and also that [2az\ + b)/d — —(2az 2 + b)/d = 
1, we can write the total particle number produced by 
the resonance for type II modes as: 



, res (II) i \ 

logn fe \z 2 ) 



1 

4^ 



Similarly, for type III modes, we obtain 

d 2 



, res (III) 
10g7J fc 



(*0 



8(-o)t 



(A6) 



(A7) 



It is easy to check that Eqs. (|A6|) and (|A7jl give the 
same result for the mode in the transition between types 
II and III, at (3 — 1, as expected. Recalling that \qo\ <C 1, 
these expressions can be approximated by those given in 

Eq. dini]). 



Hk 



(z) = -\/ az 2 +bz + c , 



(A2) 
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